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SECTION  1 


INTRODUCTION 


1.1  BACKGROUND 


Several  recent  studies  have  demonstrated  that  the  structural  effi¬ 
ciency  of  military  and  commercial  aircraft  can  be  improved  by  taking  advantage 
of  the  postbuckled  strength  of  stiffened  panels.  An  assessment  of  the  current 
postbuckled  stiffened  panel  design;  analysis  and  applications  technology 
(References  1  and  2)  shows  that  several  deficiencies  need  to  be  addressed  to 
establish  a  systematic  postbuckling  design  methodology.  In  References  1  and  3 
a  design  and  analysis  methodology  was  developed  for  flat  and  curved  stiffened 
panels  made  of  either  composite  or  metallic  materials,  and  subjected  to  either 
compression  loading  or  shear  loading.  In  practice,  however,  stiffened 
airframe  panels  are  subjecced  to  a  combination  cf  axial  compression  and  shear 
loads.  A  semi -empirical  design  methodology  for  curved  metal  panels  under 
combined  loading  exists  (Reference  4)  but  has  seen  limited  verification.  The 
present  program  was  undertaken  to  extend  the  Reference  1  and  Reference  4 
methods  for  application  to  curved  composite  panels  under  combined  uniaxial 
compression  and  shear  loading,  and  to  further  substantiate  the  metal  panel 
design  procedures. 


1.2  PROGRAM  OBJECTIVES 


The  overall  objectives  of  the  program  were  to  develop  validated 

design  procedures  and  an  analysis  capability  for  curved  metal  and  composite 

postbuckled  panels  under  combined  uniaxial  compression  and  shear  loading.  The 
c 

specific  requirements  encompassed  by  these  objectives  were  as  follows: 

*  1.  Extend  the  existing  semi -empirical  analysis  methodology  (Refer¬ 

ence  1)  into  a  design  tool  for  curved  composite  and  metal 
panels  subjected  to  combined  uniaxial  compression  and  shear 
loading.  Account  for  any  unique  failure  modes. 
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2.  Develop  a  more  rigorous  energy  method  based  analysis  to  predict 
the  displacement  and  stress  fields  in  postbuckled  panels. 

3.  Develop  a  static  and  fatigue  data  base  for  composite  and  metal 
panel  design  verification. 

4.  Develop  a  fatigue  analysis  method  for  metal  panels. 

5.  Prepare  a  procedural  design  guide.  Exercise  the  design  guide 
on  a  realistic  aircraft  component. 

The  work  performed  to  accomplish  these  objectives  is  documented  in 
this  report. 

i.3  raosRM-SJBam 


The  program  approach  and  plan  paralleled  those  in  Reference  1.  At 
the  onset,  a  technology  review  was  conducted  to  update  the  data  base  and  to 
clearly  define  the  deficiencies  in  the  static  strength,  durability  and  damage 
tolerance  design  and  analysis  of  postbuckled  metal  and  composite  panels.  The 
durability  and  damage  tolerance  technology  assessment  is  documented  in  Refer¬ 
ence  2.  As  a  result  of  this  technology  assessment,  a  semi-empirical  design 
methodology  for  curved  panels  under  combined  loading  was  established  and  a 
verification  test  program  was  planned.  In  addition,  an  energy  method  based 
approach  to  predict  the  static  response  of  postbuckled  stiffened  panels  was 
formulated.  The  technology  assessment  was  accomplished  in  Task  I  of  the 
program.  A  summary  of  the  data  gaps  identified  by  this  technology  assessment 
is  summarised  in  Table  1.  In  Task  II,  Data  Base  Development,  the  tests 
required  to  fill  the  data  gaps  identified  in  Task  I  were  conducted. 
Analytical  model  development  and  verification  was  accomplished  in  Task  III. 
Task  IV  consisted  of  technology  consolidation  where  the  results  of  the 
program  were  incorporated  in  a  Design  Guide  (Reference  5)  for  postbuckled 
structures . 

The  composite  and  metal  panels  tested  in  the  program  were  cylindric- 
ally  curved  and  identical  to  the  shear  panels  tested  in  Reference  1.  The 
design  methodology  initially  established  from  the  technology  assessment  was 
used  to  estimate  the  buckling  and  postbuckling  load  capacities  of  the  panels. 


2 


Table  1.  Summary  of  Technology  Gaps  Identified  in  Task  I. 
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Static  and  fatigue  tests  including  several  static  strain  surveys  were  conduct* 
ed  on  the  test  articles.  The  static  test  data  were  used  to  verify  the  semi- 
empirical  design  methodology,  whereas  the  fatigue  test  data  were  utilized  to 
determine  the  fatigue  failure  modes  and  obtain  load  versus  life  data  to 
formulate  fatigue  analysis  approaches.  The  test  results,  in  conjunction  with 
the  semi -empirical  design  methodology  were  used  to  update  the  Preliminary 
Design  Guide  (Reference  3) . 


1.4  REIPRT-.  OUTLINE 


This  report  details  the  correlations  between  the  results  of  the 
semi -empirical  and  rigorous  analyses,  and  the  tests  conducted.  Section  2 
describes  the  semi -empirical  analysis  methodology  and  Section  3  details  the 
development  of  the  energy  method  based  analysis.  The  actual  correlations 
between  the  analysis  and  test  results  are  presented  in  Section  4. 

Volume  I  of  the  final  report  presents  an  executive  summary  of  the 
program.  The  test  program  details  are  documented  in  Volume  II  of  the  final 
report.  Correlation  between  the  test  data  and  analyses  is  presented  in 
Volume  III.  The  Design  Guide  Update  and  the  Software  User's  Manual  are 
published  separately  as  Volumes  IV  and  V,  respectively. 


4 
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SECTION  2 


SEMI -EMPIRICAL  DESIGN  METHODOLOGY 


2.1  BACKGROUND 

An  in-depth  survey  of  the  semi -empirical  design  methods  for  metal 
panels  and  their  evolution  into  a  design  methodology  for  curved  composite 
panels  under  shear  or  uniaxial  compression  loads  is  given  in  Reference  1.  The 
shear  panel  and  compression  panel  analysis  methods  of  Reference  1  were  used  as 
the  starting  points  for  the  current  program.  Initially,  the  interaction  rules 
used  for  metal  panels  (References  4  and  6)  TTere  adopted  to  predict  buckling 
under  combined  shear  and  compression  load  g.  Test  data  were  then  used  to 
verify  these  rules  and  suggest  modifications  where  necessary.  Postbuckling 
failure  envelopes  were  developed  by  accounting  for  the  failure  modes  possible 
under  shear  loading  only,  and  under  pure  compression  loading.  Failure  predic¬ 
tions  under  combined  loading  took  into  account  load  interaction  for  stiffener 
crippling  and  skin  rupture.  This  semi -empirical  analysis  methodology  is  de¬ 
tailed  in  the  following  subsections. 


2.2  DESIGN  METHODOLOGY 


A  complete  static  analysis  of  postbuckled  structures  consists  of 
predicting  the  initial  buckling  loads,  the  failure  or  ultimate  load  of  the 
panel  after  buckling,  and  the  local  skin  and  stiffener  displacement  and  stress 
fields.  The  latter  predictions  are  especially  required  for  metal  panel 
fatigue  analysis.  The  semi-empirical  methodology  detailed  in  this  section  can 
be  used  to  obtain  the  initial  buckling  and  the  failure  loads.  The  energy 
method  based  analysis  described  in  Section  3  is  useful  in  predicting  the  local 
stresses  and  displacements. 

The  semi -empirical  analysis  method  was  selected  as  a  design  tool  for 
postbuckled  structures  to  provide  a  quick,  inexpensive,  and  reasonably  accu¬ 
rate  but  conservative  design  methodology.  The  scope  of  this  program  encom- 
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passed  cylindrically  curved  stiffened  panels  loaded  in  uniaxial  longitudinal 
compression  and  shear.  Since  the  Reference  1  methodology  is  the  basis  for  the 
combined  loading  design  procedure,  the  semi -empirical  analysis  described  in 
this  section  applies  to  cylindrically  curved  panels  under  simultaneously 
acting  longitudinal  compression  and  in-plane  shear. 

The  essence  of  the  combined  loading  design  procedure  is  summarized  „ 

in  Figure  1.  As  can  be  seen  in  the  figure,  the  curved  panel  is  analyzed  for 
compression  and  shear  loads  independently  according  to  Reference  1  methods.  * 

Buckling  loads  under  combined  loading  are  predicted  using  the  parabolic 
interaction  rule  developed  for  metal  panels  (References  4  and  6) .  Failure 
analysis  requires  consideration  of  failure  modes  under  shear  and  compression 
acting  independently  and  those  due  to  the  interaction  of  the  loads.  The 
failure  modes  affected  by  combined  loading  are  stiffener  crippling,  and  skin 
rupture  under  tensile  loading  determined  from  a  principal  strain  analysis  and 
the  maximum  strain  criterion.  The  following  paragraphs  present  the  detailed 
equations  necessary  for  a  semi -empirical  analysis. 


2.2.1  Skin  Bu-9kUng..S.train/!koad 


The  shear  and  compression  buckling  strains  or  loads  for  metal  and 
composite  panels  are  calculated  according  to  the  equations  given  in  Reference 
1.  For  continuity  these  equations  are  summarized  below. 

The  compression  buckling  stress  for  curved  metal  sheet  panels  can 
be  calculated  from: 


where , 


fCR 


K  jr2E 
c _ 

12(l-i/2) 


2 


For  buckling  stress,  psi 

thickness  of  the  skin,  in 

bs  stiffener  spacing  measured  between  fastener  lines,  in 


(1) 


6 


Figure  1.  Semi-empirical  Analysis  Approach  for  Postbuckled 
Stiffened  Panels  Under  Combined  Loading. 
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E ,v  modulus  and  Poisson's  ratio  for  the  sheet  material 

Kc  buckling  coefficient  determined  from  Figure  2  (References  6  and 

7) 

The  theoretical  value  of  Kc  is  obtained  from  the  buckling  equations  for  thin 
cylindrical  shells  and  is  a  function  of  the  nondimens ional  curvature  Z  of  the 
panel  expressed  as 

b2  (l-t/2)h 
Z  -  -S - 

rtw 

where  r  is  the  radius  of  the  cylindrical  panel.  Experimental  data  (Reference 
7)  have  shown  that  Kc  is  also  a  function  of  the  r/t  ratio  for  the  panel.  The 
design  curves  of  Figure  2,  obtained  from  test  data,  show  this  dependence  of  Kc 
on  r/t. 

Compression  buckling  strains  for  curved  composite  panels  can  be 
accurately  determined  through  the  use  of  computer  codes  SS8  (Reference  8)  and 
BUCLASP-2  (Reference  9),  for  example.  However,  for  an  approximate  calculation 
of  the  skin  buckling  strain,  the  simplified  equation  given  below  can  be  used. 


*  ' 

2 

* 

'  ' 

2 

9  y 

nwr 

1 

Du  +  2(d12+2d66) 

nL 

+  D22 

nL 

L 

. 

^xw^ 

mbw 
L  WJ 

mbw 
L  WJ 

+ 


E 


m* 

L 


R2 


®x w 


E  E 

’  i 

2 

'  ' 

O..  1?  xw  vw 

^xyw^yw  ~ — 

bxyw 

nL 

mbw 

+  Eyw 

nL 

mbw 

(2) 


where  Djj  are  the  terms  of  the  bending  stiffness  matrix  of  the  composite  skin, 
Exw,  Eyw ,  Gxyw>  uxyw>  atu*  tw  are  the  web  elastic  constants  and  thickness, 
respectively,  L  is  the  panel  length,  bw  is  the  width  of  the  skin,  R  is  the 
radius  of  curvature  of  the  panel  and  n  and  m  are  integer  coefficients  repre¬ 
senting  the  number  of  half  buckle  waves  in  the  width  and  length  direction, 
respectively.  The  lowest  value  of  strain  for  various  values  of  n  and  m 
represents  the  buckling  strain  of  the  specimen. 
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Axial  Compression  Buckling  Coefficients  for  Long  Curved  Plates 
(Reference  6). 


The  effective  width  of  the  skin,  bw,  was  assumed  to  be  equal  to  the 
distance  between  the  two  adjacent  stiffeners  measured  from  one  stiffener 
flange  centroid  to  the  next  stiffener  flange  centroid.  Note  that  bw  is  less 
than  the  stringer  spacing  hs. 

Equation  2  was  derived  in  Reference  10  from  the  equations  developed 
for  the  buckling  of  orthotropic  complete  cylinders  by  making  simplifying 
assumptions . 


The  shear  buckling  stress  or  strain  for  composite  webs  can  be  calculated  using 
program  SS8  (Reference  9) .  The  buckling  stress  for  curved  metal  webs  can  be 
calculated  usinj  * 


fcr, elastic 


K  .n2  Eh2 
si _ s 

12R2Z2 


12R2Z2 


if  hr  >  hs 


if  hs  £  hr 


(3) 


where , 


Ksl>Ks2 


R 

E 

Z 


v 


critical  shear  stress  coefficients  for  simply  sup¬ 
ported  curved  plates,  given  in  Reference  6 

panel  radius,  in. 

Young's  modulus  for  the  material,  psi 


h2 

Rtw 


Ja-uh 


if  hr  >  hs 


h2 


Rt,, 


if  hs  >  hr 


Poisson's  ratio  for  the  material 


The  composite  panel  buckling  loads  obtained  from  program  SS8  are  in 
terns  of  running  loads  Nxcr  and  NXyCr  for  compression  and  shear  loading, 
respectively.  The  critical  buckling  stresses  for  metal  panels,  Equations  1 
and  3  can  be  converted  to  running  loads  as  follows: 
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o 

^xcr  ■*  ^cr '  t-w 
o 

^xycr  “  Tcr,  elastic  ‘  t-w 
where  tw  is  the  skin  thickness. 

For  combined  compression  and  shear  loading,  the  buckling  loads  can 
be  computed  from  (ReA.rence  6): 

2 

Rc  +  Rs  -  1  (4) 

o  o  o  o 

where,  Rc  -  ^xcr/^xcr  and  Rs  ^xycr/Nxycr .  ^xcr  and  ^xycr  are  the  pure 

compression  and  pure  shear  buckling  loads,  respectively,  and  ®xcr  and  ^xycr 

are  the  buckling  loads  when  the  shear  and  compression  loads  are  acting 

simultaneously.  The  presence  of  compression  stresses  reduces  the  shear 

buckling  stress  and  vice  versa. 


2.2.2  Failure  Analysis  and  Margin  Computation 


Failure  analysis  of  postbuckled  structures  requires  identification 
of  all  possible  failure  modes  and  calculating  the  loads  corresponding  to  the 
critical  failure  mode.  For  curved  panels  under  combined  loading  a  failure 
envelope  spanning  the  load  ratio  Nx/NXy  values  of  0  (i.e.,  Nx-0,  NXy*0)  to  <* 
(i.e.,  Nx*0,  NXy-0)  is  a  convenient  means  for  identifying  the  critical  failure 
mode.  The  procedure  to  develop  this  failure  envelope  is  detailed  in  the 
following  subsections. 


2 . 2 . 2 . 1  Compression  Loading  Failure  Analysis  (Nx/NXy-«>.  Nyy-0) 

The  analysis  for  failure  under  compression  loading  has  been  devel¬ 
oped  and  documented  in  Reference  1.  Under  compression  loading  the  possible 
failure  modes  are: 


11 


1.  Euler  buckling  of  the  stiffened  panel 

2.  Stiffener  crippling 

3.  Stiffener/ skin  separation  for  composite  panels  with  cocured  or 
bonded  stiffeners 

4.  Skin  permanent  set  for  metal  panels. 

Euler  Buckling  Strain  Calculations .  The  Euler  buckling  strain  for  a 
stiffened  panel  is  calculated  by  treating  the  panel  as  a  wide  column  with  the 
width  set  equal  to  the  stiffener  spacing.  The  critical  strain  is  calculated 
using  the  standard  column  equation: 


€CR 


2 

C*  El 
EA  L2 


(5) 


where,  El  is  the  equivalent  bending  stiffness  of  the  panel,  EA  is  the  equiva¬ 
lent  axial  stiffness,  L  is  the  panel  length,  and  C  is  the  end  fixity  coeffi¬ 
cient.  The  fixity  coefficient  depends  upon  the  support  conditions  at  the 
panel  ends.  Most  compression  panels  are  tested  by  flat  end  testing  and  the 
results  obtained  by  using  C  -  4  are  unconservative;  therefore,  a  value  of  C  - 
3  is  recommended.  The  values  of  C  for  other  end  conditions  can  be  obtained 
from  Reference  6  (Subsection  A18.23). 

Stiffener  Crippling  Strain/Stress  Calculation.  The  crippling 
strength  of  metal  stiffeners  is  calculated  using  the  well  established  Needham 
or  Gerard  methods  documented  in  Reference  7.  In  the  present  program,  the 
Gerard  method  was  used  since  it  is  a  generalization  of  the  Needham  method  and 
was  derived  from  a  broader  data  base.  The  empirical  Gerard  equation  for 
calculating  the  crippling  stress  for  2  corner  sections,  such  as  the  Z,  J  and 
channel  sections ,  is : 


F 

_cs 

Fcy 


3.2 


■ 

2 

'  ’ 

1/3- 

tf 

E 

A 

»  . 

.Fcy. 

0.75 


(6) 


where 


Fes  “  crippling  stress  for  the  section,  psi 

FCy  -  compressive  yield  stress  of  the  material,  psi 

t  -  element  thickness,  in. 

A  -  section  area,  in^ 

A  design  curve  based  on  Equation  6  is  shown  in  Figure  3  taken  from  Reference 
6  Additional  crippling  equations  that  apply  to  sections  other  than  2  corner 
sections  are  also  given  in  Reference  6. 

In  order  to  calculate  the  crippling  strains  for  stiffeners  made  of 
composite  materials,  a  semi-empirical  methodology  was  developed  in  the  pro¬ 
gram.  The  methodology  consists  of  modeling  the  stiffener  in  terms  of  inter¬ 
connected  flat  plate  elements,  calculating  the  initial  buckling  and  crippling 
strains  for  each  element,  and  determining  the  crippling  strain  for  the  stiff¬ 
ener  as  the  lowest  strain  that  causes  crippling  of  the  most  critical  element 
in  the  stiffener  section.  It  should  be  noted  here  that  the  absolute  minimum 
of  the  crippling  strains  for  the  various  plate  elements  is  not  necessarily  the 
stiffener  crippling  strain;  element  criticality  with  respect  to  stiffener 
stability  has  to  be  considered  as  well.  The  procedural  details  of  this 
methodology  given  in  the  following  paragraphs  provide  additional  clarifica¬ 
tions  relating  to  the  determination  of  the  most  critica1  plate  element. 

The  first  step  in  calculating  the  stiffener  crippling  strain  is  to 
model  the  stiffener  as  an  inuerconnected  assembly  of  plate  elements.  As 
examples,  plate  element  models  of  a  hat-section  and  a  J-section  stiffener  are 
shown  in  Figure  '+.  The  hat-seccion  stiffener  is  made  up  of  four  elements, 
whereas,  the  J-section  stiffener  consists  of  nine  elements. 


The  crippling  strains  for  the  plate  elements  are  calculated  from 
empirical  equations  of  the  form 


(7) 


where 
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.1  1.5  2  2.5  3  4  5  6  7  8  8  10  15  20  25  30  40  50  6070  8090100 

Figure  3.  Crippling  Stress  Fcs  for  Two  Corner  Sections  e.g.,  Z,  J, 
and  Channel  Sections  (Reference  6,  Figure  C7-9), 
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Figure  4.  Plate  Element  Models  of  Hat-  and  J-Section  Stiffeners. 
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ecs  «  crippling  strain  of  the  plate  element 

ecr  -  initial  buckling  strain  of  the  plate  element 

ecu  “  compression  ultimate  strain  for  the  plate  element  laminate 

a,fi  -  material  dependent  coefficients  obtained  from  test  data 


Equation  7  has  the  same  functional  form  as  that  used  by  Gerard  (Reference  7) 
for  metal  stiffeners.  The  coefficients  a  and  /3  depend  on  the  plate  edge 
conditions  and  have  been  obtained  in  References  11  and  12  from  a  large  data 
base  for  plate  elements  that  are  connected  on  both  sides  (e.g.,  Elements  2,  3, 
and  4  of  the  hat-section  stiffener  shown  in  Figure  4).  The  crippling  strain 
for  stiffener  plate  elements  connected  on  both  sides  is  given  by  (Reference 
12): 


€cs 


0.56867ecr 


0.47567 


(8) 


where  ecr,  the  buckling  strain  for  the  plate  element,  is  given  by  (Reference 
13): 


ecr 

In  Equation  9 
b 
t 

Ex 


2n 


b2t  Ex 


a/°11d22  +  d12  +  2d66 


(9) 


plate  element  width 
plate  element  thickness 

compression  modulus  of  the  plate  laminate  along  the 
longitudinal  direction 

terms  from  the  laminate  bending  stiffness  matrix,  (i,j  - 

1,  2,  6) 


Equation  9  applies  to  plate  elements  for  which  the  length- to-width  ratio  (L/b, 
where  L  -  stiffener  length)  is  at  least  4. 


The  crippling  strain  for  plate  elements  that  are  connected  on  one 
side  only  is  calculated  using  the  following  equation: 


■cc 


0.4498c 


cr 


10.72715 


cu 


-cr 


(10) 
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where , 


«cr 


b^t  L^t  Ex 


L  *•  length  of  the  stiffener 


(ID 


with  the  other  nomenclature  remaining  the  same  as  for  Equations  8  and  9 . 

The  coefficients  in  Equation  10  were  obtained  by  fitting  Equation  7 
to  the  crippling  data  generated  from  tests  on  one-edge  free  plates  in  Refer¬ 
ences  11  and  12.  Data  for  two  material  systems,  T300/5208  and  AS/3501  graph¬ 
ite/epoxy,  were  pooled  to  obtain  Equation  10. 


In  Equations  8  through  11,  the  thickness  of  plate  elements  attached 
to  the  skin  is  taken  as  the  sum  of  the  plate  element  and  the  cocured  skin 
thicknesses.  In  the  case  of  the  hat- section  stiffener,  crippling  strains  for 
plate  elements  representing  the  skin  only,  such  as  Element  4  in  Figure  4  are 
also  calculated.  Another  consideration  in  calculating  the  crippling  strain 
for  stiffener  flange  elements  attached  to  the  skin  is  the  choice  of  an  appro¬ 
priate  element  width.  For  example,  in  most  practical  designs  the  stiffener 
flanges  attached  to  the  skin  are  tapered  by  dropping-off  plies  as  shown  in 
Figure  5  for  a  hat-section  stiffener.  The  flange  plate  element  width  in  this 
case  Is  defined  as  the  width  to  the  end  of  the  taper  with  the  weighted  average 
of  the  element  thickness  added  on  to  the  attached  skin  thickness  to  obtain  the 
total  thickness  for  use  in  Equations  8  through  11. 


Equations  8  through  11  are  quite  general  in  nature  and  take  into 
account  ply  composition,  stacking  sequence,  and  material  characteristics.  The 
ply  composition,  i.e.,  the  percentages  of  0°,  45°,  and  90°  plies,  is  reflected 
in  the  compression  ultimate  strain  ecu.  Stacking  sequence  effects  are  ac¬ 
counted  for  in  the  expression  for  ecr  where  the  bending  stiffnesses  Dj_j  are 
used.  The  Djj's  an<*  ecu  also  account  for  mechanical  property  changes  from  one 
material  system  to  another.  Use  of  strair  rather  than  stress  for  crippling 
calculations  provides  another  significant  advantage  in  that  laminate  non¬ 
linearity  (e.g.,  stress-strain  response  of  ±45"  laminates)  is  accounted  for  by 
way  of  the  compression  ultimate  strain  ecu. 
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Figure  5.  Ply  Drop-Offs  in  Hat-Section  Stiffener. 
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Failure  Load  Calculation.  The  failure  load  for  the  panel  is  deter¬ 
mined  as  the  lowest  of  the  loads  calculated  for  the  various  instability  modes 
mentioned  above,  for  stiffener-web  separation  in  composite  panels,  and  for 
skin  or  stiffener  yielding  in  metal  panels.  The  methods  for  failure  load 
calculation  are  given  in  the  following  paragraphs. 

Failure  Load  Due  to  Euler  Buckling.  The  failure  load  due  to  Euler 
buckling  is  calculated  using  the  following  equation: 

E 

Pe  “  €cr  (^xs^s  +  ^xwbwtw)  (12) 

where , 

E 

ecr  -  Euler  buckling  strain  determined  using  Equation  5 

Ejjg  -  Compression  modulus  of  the  stiffener  in  the  loading 

direction 

As  -  Cross-sectional  area  of  the  stiffener 

Ejh,  -  Compression  modulus  of  the  web  (skin)  in  the  loading 

direction 

bw  -  Stiffener  spacing 

tw  -  Skin  thickness 

Failure  Load  Due  to  Stiffener  Crippling.  In  order  to  determine  the 
failure  load  due  to  stiffener  crippling,  it  is  necessary  to  determine  the  load 
carried  by  the  stiffener  and  the  panel  web  individually.  The  load  carried  by 
the  stiffener  (Ps)  is  determined  as  follows: 

1.  Determine  the  two  lowest  crippling  strains  (eCcl)  and  (ecc2)  °f 
all  the  elements  making  up  the  cross-section  using  Equations  8 
through  11. 

2.  If  the  element  with  the  lowest  crippling  strain  ( eccl ) 
normal  to  the  axis  of  least  bending  stiffness  of  the  cross- 
section,  the  stiffener  will  fail  at  a  strain  equal  to  ecci,  and 
the  corresponding  failure  strain  of  the  stiffener  is  given  by: 

PS  ”  ^xs^x  eccl 
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3.  If  the  element  with  the  lowest  crippling  strain  is  parallel  to 
the  axis  of  least  bending  stiffness  of  the  cross-section,  the 
stiffener  will  carry  additional  load  until  the  second  member  in 
the  cross-section  becomes  critical  due  to  crippling.  In  this 
case  the  load  carried  by  the  stiffener  is  given  by: 

Ps  —  (EA)]_  (eCcl  “  €cc2)  +  ecc2  ®xs^s  (1^) 

where  (EA)i  is  the  extensional  stiffness  of  the  member  becoming 
critical  first,  and  the  stiffener  failure  strain 


cc 

€s  ”  €cc2 

The  total  load  carried  by  the  panel  is  the  sum  of  the  load  carried 
by  the  stiffener  up  to  crippling  and  the  load  carried  by  the  buckled  skin.  In 
order  to  calculate  the  load  carried  by  the  skin,  the  effective  width  concept 
is  utilized.  The  effective  width  for  metal  panels  is  calculated  using  the 
semi-empirical  equation  given  below  (Reference  6) : 


w 


where 


(15) 


w  -  effective  width  of  the  skin  after  initial  buckling 
-  skin  thickness 
Fst  ~  stress  in  the  stringer 

For  composite  panels,  in  the  absence  of  any  other  guidelines, 
Equation  15  expressed  in  terms  of  strain  is  used  to  compute  the  effective  skin 
width .  Thus , 


-0.5 

w  -  1.9tw(es)  (15A) 

for  composite  skins  where  es  -  strain  in  the  stiffener. 

Thus,  the  total  load  carried  by  the  panel  for  a  stiffener  crippling 
mode  of  failure  is  given  by: 

PCc  *“  Ps  +  Pw  (16) 
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where 

Pcc  -  load  carried  by  the  panel  at  stiffener  crippling 

Ps  -  stiffener  load  given  by  Equation  14 

Pw  "  load  carried  by  the  skin 
The  load  Pw  is  calculated  as: 


Pw  *■  Fcs  w  tw  —  1.9tw  yEFcs  (17) 

for  metal  panels,  and  for  composite  panels  as: 

2  cc  0-5 

Pw  *•  1.91^,  Ej-W( e s  )  (18) 

Failure  Load  Due  to  Stiffener /Web  Separation.  Failure  of  composite 
stiffened  panels  due  to  stiffener/web  separation  is  a  common  mode  of  failure 
in  the  postbuckling  range.  It  is  extremely  difficult  to  predict  this  failure, 
even  by  using  rather  sophisticated  analysis  methods.  The  attempts  to  date  on 
making  such  predictions  have  been  inconclusive.  A  simple  empirical  equation 
to  predict  such  failure  was  developed  in  this  program.  The  correlation  of 
experimental  data  with  the  predicted  failure  loads  based  upon  this  equation  is 
surprisingly  good.  The  empirical  equation  was  derived  by  analogy  with  the 
crippling  data  for  plates  with  one  edge  simply  supported  and  one  edge  free. 
It  is  hypothesized  that  when  the  panel  web  strain  reaches  the  crippling  strain 
the  interfacial  stresses  become  high  enough  to  cause  failure.  The  equation 
should  represent  the  lower  bound  on  predicted  failure  loads.  Any  attempts  to 
improve  the  interface  (for  example,  by  stitching,  riveting,  etc.)  can  result 
in  higher  failure  loads. 

?ss  "  €ss  (Exs^S  +  Fxwbwtw)  (19) 


where 


0.4498  ecr  -SB 


0.72715 


€ ss  “  Failure  strain  for  stiffener/web  separation 

Pss  ”  Failure  load  for  the  stiffener/web  separation  mode 
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The  metal  compression  panel  analysis  methodology  outlined  in  the 
preceding  paragraphs  has  been  experimentally  validated  (e.g.,  Reference  7)  and 
is  representative  of  current  usage.  In  the  case  of  composite  panels, 
available  composite  compression  panel  test  data  were  utilized  to  validate  the 
semi -empirical  analysis  (Reference  1) . 


V 

2 . 2 . 2 . 2  Shear  Loading  Failure  Analysis  (Nx/NKy-0:  N^y*0 ) 

* 


t 
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Flat  or  curved  shear  panel  analysis  is  accomplished  by  means  of  the 
semi -empirical  tension  field  theory  developed  by  Kuhn  (Reference  4)  for  metal 
panels.  In  Reference  1  the  tension  field  theory  was  modified  for  application 
to  composite  shear  panels  by  taking  into  account  material  anisotropy. 

The  essential  elements  of  the  generalized  (for  application  to  metals 
as  well  as  composites)  tension  field  theory  and  its  application  are  summarized 
in  Figure  6.  Details  of  the  semi -empirical  analyses  required  to  perform  the 
various  steps  in  Figure  6  are  given  in  the  following  paragraphs.  The  equa¬ 
tions  as  presented  below  pertain  to  cylindrically  curved  composite  panels  and 
to  flat  composite  panels  if  terms  incorporating  the  radius  of  curvature  R  are 
set  equal  to  zero.  Use  of  the  appropriate  values  for  elastic  constants  in  the 
equations  permits  their  direct  application  to  metal  panels.  The  analysis 
procedure  is  based  entirely  on  the  theory  presented  in  Reference  4  unless 
specifically  noted. 


Computation  of  the  Diagonal  Tension  Factor.  The  diagonal  tension 
factor  k  characterizes  the  degree  to  which  diagonal  tension  is  developed  in 
the  skin  of  stiffened  panels  loaded  in  shear.  A  value  of  k  -  0  characterizes 
an  unbuckled  skin  with  no  diagonal  tension;  a  value  of  k  -  1.0  characterizes  a 
web  in  pure  diagonal  tension.  The  diagonal  tension  factor  is  computed  using 
the  following  expression: 
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COMPUTE  DIAGONAL  TENSION 
FACTOR  ’K' 
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Figure  6.  Application  of  Tension  Field  Theory  to  Shear  Panels. 
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where 


i 


i. 


K 


ty  -  web  thickness 

hr  -  ring  spacing 

hs  -  stringer  spacing 

R  -  panel  radius 

t  -  applied  shear  stress  -  Nxy/t^ 

o 

rcr  -  buckling  shear  stress  of  web  under  pure  shear  conditions 

NxycrAw 

The  pure  shear  buckling  stresses  for  composite  and  metal  panels  are  calculated 
using  the  techniques  given  in  Subsection  2.2.1. 


Computation  of  Diagonal  Tension  Angle  *a* .  An  initial  value  is 
assigned  to  the  diagonal  tension  angle  'a'  that  defines  the  angle  of  the 
'folds'  in  the  buckled  skin.  For  curved  web  systems  a-30°  was  found  to  be  a 
convenient  starting  point.  The  actual  value  of  a  is  determined  by  the  itera¬ 
tive  procedure  outlined  below. 


Using  the  assumed  initial  value  of  a  a  'new'  value  for  a  is  calcu¬ 
lated  by  the  equation: 
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if  hr  >  hs 


Rf 


lRJ 


1 

8 


2 

Tan^a 


For  eccentric  stringers  and  rings 


if  hs  >  hr 


(22d) 


EAr 


(22e) 


In  Equations  22,  e  is  the  skin  strain  in  the  diagonal  tension  direction,  and 


c g  and  «r  are  the  strains  in  the  stringer  and  the  ring  leg  attached  to  the  web 
averaged  over  their  lengths,  respectively.  E^*,  Evs,  and  Ewr  are  the  web 
moduli  in  the  direction  of  the  tension  field,  stringers  and  rings,  respec¬ 
tively.  Grs  is  the  web  shear  modulus,  EAS  and  EAr  are  the  effective  axial 
stiffnesses  of  the  stringers  and  the  rings,  respectively,  calculated  with 
respect  to  the  skin  mid-surface.  El  is  the  bending  stiffness  about  the  stiff¬ 
ener  neutral  axis  and  El  the  bending  stiffness  about  the  web  midsurface. 

In  general,  oj_,  the  new  diagonal  tension  angle  will  not  equal  the 
initially  assumed  value  of  30°.  Therefore,  is  used  as  the  next  guess  and 
the  computations  of  Equation  22  are  repeated  until  the  process  converges, 
i.e. ,  anew  =  «old’ 

Once  the  diagonal  tension  angle  has  been  determined  with  sufficient 
accuracy,  the  next  step  is  to  compute  the  margins  of  safety. 

Computation  of  Stringer  and  Frame  Margins  of  Safety .  The  diagonal 
tension  angle  value  computed  above  is  now  substituted  in  Equations  22  to 
obtain  the  diagonal  tension  strain  in  the  skin,  the  stringer  strain,  and  the 
ring  strain.  Next,  the  stringer  and  ring  strains  averaged  over  the  cross 
section  and  the  length  (eave)  and  the  maximum  strains  in  the  legs  attached  to 
the  web  (6max)  ate  computed  using  the  following  equations: 
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The  stringer  and  ring  crippling  mode  of  failure  is  then  analyzed  for  by 
computing  the  stringer  and  ring  forced  crippling  strains  (eos  and  €or,  respec¬ 
tively)  using  the  following  equations: 
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where  ean  is  the  laminate  allowable  strain,  Ecs  and  Ecr,  are  the  compression 
modulus  of  the  stringer  and  ring  leg  attached  to  the  web,  respectively,  and 
tus  and  tur  are  the  thickness  of  the  stringer  and  the  ring  leg  attached  to  the 
web. 


The  critical  stiffener  strains  corresponding  to  the  bending  stiff¬ 
ness  required  for  stiffener  stability  are  calculated  using  Equations  29  and 
30. 
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(29) 
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where  €sg  and  rg  are  the  Euler  buckling  strains  for  the  stiffener  and  the 
ring,  respectively. 

The  margins  of  safety  can  now  be  computed  for  each  of  the  possible 
failure  modes  by  comparing  the  calculated  strain  values  with  the  allowables. 
Thus,  to  ensure  positive  margins,  the  following  failure  modes  are  examined  and 
the  corresponding  inequalities  verified. 


1.  For  stringer  and  ring  stability 
i . e . ,  no  column  failure 


£sB  >  £sj 
£rB  >  6rc 


2.  For  stability  of  the  entire  EIS  >  Estw( — §•  -  2)h£ 


panel,  i.e.,  to  prevent  buck¬ 
ling  of  the  Web  as  a  whole,  be¬ 
fore  formation  of  the  tension 
field 


3  3h 

EIr  >  Ertw( — £  -  2)hr 


3.  For  prevention  of  forced  cripp¬ 
ling  of  stiffeners 


£os  >  £smax 
£or  >  £rmax 


An  additional  check  needs  to  be  performed  for  metal  panels  where 
yielding  or  permanent  set  in  the  web  is  likely  due  to  excessive  skin  deforma¬ 
tion.  The  only  available  criterion  for  permanent  set  check  has  been  empiri¬ 
cally  obtained  from  tests  on  flat  aluminum  metal  panels.  Its  applicability  to 
other  materials  or  curved  panels  has  not  been  verified.  Thus,  in  the  absence 
of  any  other  guidelines,  the  flat  panel  requirement  that  the  maximum  allowable 
value  of  the  diagonal  tension  factor  ka^  be  limited  to 

kall  “  °-?8  -  (t-0.012)0-50  (32) 

at  design  ultimate  load  to  prevent  permanent  buckling  of  the  web  at  limit 
load,  is  used  in  the  present  analysis. 
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2 . 2 . 2 . 3  Combined  ..Loading  Failure  Analysis  (Nv/Nyy-B:  Nrcr^Nyycr-A) 


The  effects  of  shear  and  compression  loading  interaction  have  to  be 
accounted  for  in  a  combined  loading  failure  analysis.  For  the  combined 

loading  case,  the  additional  considerations  are: 

1.  The  buckling  stresses  are  reduced  in  accordance  with  the 
interaction  given  in  Equation  4. 

2.  Compression  stresses  in  the  stiffeners  prior  to  buckling  are 
those  due  to  the  directly  applied  compression  only.  However, 
after  buckling  the  compression  stresses  due  to  diagonal  tension 
must  be  added  to  the  direct  compression. 

3.  The  allowable  stress  calculation  for  the  stiffeners  must 

account  for  an  interaction  between  the  forced  crippling  (panel 
shear  induced)  and  natural  crippling  (direct  compression 

induced)  modes  of  stiffener  failure. 

4.  Calculation  of  the  stiffener  stresses  due  to  applied  shear 

loads  is  modified  to  account  for  the  presence  of  the  compres¬ 
sion  load. 

The  buckling  interaction  equation  can  be  rewritten  as 

°  °  9 

Nxcr/^xcr  +  (NXycr/NXyCr)  -  1 

then, 
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^xycr  ”  NXycr  J 1- (Nxcr/Nxcr) 

The  diagonal  tension  factor  k  is  expressed  as 
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where 


NXy  ”  applied  shear  load 

NXyCr  “  shear  buckling  load  for  combined  loading  as  calculated 
from  Equation  33. 

Calculation  of  k  using  Equation  34  is  subject  to  the  auxiliary  rules  that  if 
hs  >  hr  then  replace  hr/hs  with  hs/hr  and  if  the  resulting  ratio  is  greater 
than  2  then  use  a  value  of  2  for  the  ratio. 

The  diagonal  tension  angle  a  is  computed  iteratively  using  the  same 
procedure  as  for  pure  shear,  but  with  appropriate  modifications  to  the 
stiffener  strain  expressions.  Thus,  if  a-30°  initially  then  the  new  a  is 
calculated  using  Equation  22  where  e ,  the  skin  strain  is  obtained  from 
Equation  22a,  the  ring  strain  from  Equation  22c  and  the  stiffener  strain  from 
the  following  expression: 


£g  -  - ikr-Cptot -  (35) 

Tea 

r-S-  +  o.sa-kjEwgRg 

hgtw 

2 

where  Rc+Rs  -  1. 

As  in  the  pure  shear  case,  sufficient  iterations  need  to  be  per¬ 
formed  so  that  crnew  53  aQid. 

Computation  of  Stiffener  Margin  of  Safety.  The  total  stiffener  load 
can  be  expressed  as: 

Ps  “  Px  +  Pxy  (36) 

where  Px  is  the  load  in  the  stiffener  due  to  direct  compression  and  PXy  is  the 
load  in  the  stiffener  due  to  the  diagonal  tension  folds.  The  resulting 
stiffener  strain  can  be  expressed  as  (References  6  and  7) 
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where  the  negative  signs  denote  a  compression  strain,  w  is  the  effective  width 
of  the  skin  after  buckling  as  obtained  from  Equations  17  or  18  and  Rs  from 
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Equation  4.  The  average  and  maximum  strains  in  the  stiffener  can  be  computed 
by  analogy  to  the  pure  shear  case ,  i . e . , 
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In  computing  margins  of  safety  for  stiffener  design,  the  above 
strains  have  to  be  compared  against  the  Euler  buckling  strain  and  the  stiffen¬ 
er  crippling  allowable  strain.  For  Euler  buckling,  it  is  immaterial  whether 
the  stiffener  compressive  strain  arises  from  the  direct  compression  load,  Px, 
or  from  the  diagonal  tension  action  caused  by  PXy.  Ec  1  <tx  buckling  failure  is 

assumed  to  take  place  when  «s  given  by  Equation  38  above  reaches  egg  given 

ave 

by  Equation  29.  The  nature  of  stiffener  crippling  under  combined  loading, 
however,  requires  that  the  interaction  between  the  strain  due  to  direct  com¬ 
pression  and  the  strain  due  to  diagonal  tension  be  accounted  for.  This  is 
because  crippling  under  diagonal  tension  is  caused  by  forced  deformation  of 
the  stiffener  leg  attached  to  the  web,  whereas  direct  compression  causes  crip¬ 
pling  failure  by  local  instability  of  the  entire  stiffener  section.  An  empir¬ 
ical  expression  for  this  interaction  has  been  given  in  Reference  4  for  curved 
metal  panels.  For  generic  application  to  metal  and  composite  panels  the 
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Reference  4  interaction  is  expressed  in  terms  of  strains  as  follows: 
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where 


co 

es  -  the  direct  compression  strain 


so 

€ g  -  the  compression  strain  due  to  diagonal  tension  which  cause 
stiffener  crippling  v;hile  acting  simultaneously 


cc 

«s  -  the  stiffener  crippling  strain  under  pure  compression 

loading  as  computed  from  Equations  9  and  10,  and  the 
procedure  given  in  Section  2.2.2. 

€os  “  the  forc®d  crippling  strain  of  the  stiffener  under  pure 

shear  loading  calculated  from  Equation  27. 


The  margin  of  safety  is  computed  as  follows: 
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Computation  of  Ring  Margin  of  Safety.  Metal  panel  test  data  show 
that  the  hoop  compression  stresses  in  the  ring  due  to  diagonal  tension  are 
unaffected  by  the  axial  compression  on  the  curved  panel  as  a  whole.  There¬ 
fore,  Equations  22c,  25,  and  26  can  be  readily  used  to  compute  the  ring 
strains  and  margins  of  safety. 
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2.2.3  Automated  Semi -Empirical  Design  Methodology 


The  design  procedure  outlined  above  has  been  coded  in  a  computer 
program  called  PBUKL  for  use  as  a  design  tool.  Detailed  documentation  of  this 
program  is  given  in  Reference  14.  The  program  is  an  extension  of  TENWEB, 
works  interactively,  and  has  several  built-in  stiffener  profiles  for  design 
flexibility. 

Program  PBUKL  was  used  to  design  the  curved  panels  tested  in  this 

s  tudy , 
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SECTION  3 


ENERGY  METHOD  BASED  ANALYSIS  DEVELOPMENT 

3.1  PROBLEM  FORMULATION 

The  energy  approach  was  used  for  problem  formulation.  The  problem 
was  formulated  for  a  stiffened,  curved  anisotropic  laminated  plate.  The 
laminate  was  assumed  to  be  balanced  and  symmetric.  A  small  imperfection  in 
the  lateral  displacement  was  also  included  in  the  formulation.  The  panel 
geometry  and  the  coordinate  system  are  shown  in  Figure  7.  This  figure  also 
shows  the  relationship  between  the  overall  postbuckling  structural  configura¬ 
tion  and  the  panel  geometry  used  in  the  analysis.  Since  the  adjacent  bays  are 
assumed  to  deform  in  an  identical  fashion,  a  single  bay  was  analyzed.  Figure  7 

•k  k  * 

shows  that  the  material  properties  for  the  skin  are  Af  j  and  Djj  ,  where  Af  j 
*  *  *  *  *  *  *  * 
(A^,  A^2>  ^22  an<i  Agg)  is  the  skin  stiffness  matrix  and  Dfj  (D^,  Df2>  °16> 

it  it  A 

^22*  ®26  an<*  Dgg)  *-s  t*ie  skin  rigidity  matrix.  The  material  axes  1  and  2  are 
assumed  to  coincide  with  the  panel  geometry  coordinate  axes  x  and  y, 
respectively.  The  panel,  with  length  a  and  width  b,  is  bounded  by  stringers 
along  the  straight  edges  and  frames  or  rings  along  the  curved  edges .  The 
cross-sectional  area.  Young's  modulus  and  moment  of  inertia  of  the  stringers 
are  As>  Es  and  Is  respectively.  Those  of  the  frames  are  Af,  Ef  and  If.  The 

radius  of  the  curved  panel  is  R. 

The  energy  expressions  are  written  in  terms  of  the  displacement 
components  u,  v  and  w  in  the  x,  y  and  z  directions,  respectively.  The  panel 
is  assumed  fixed  along  x  -  o  and  subjected  to  a  system  of  combined  compression 
(Nx)  and  shear  (NXy)  load  along  the  edge  x  -  a.  The  boundary  conditions  for 
the  displacement  components  are  therefore  given  by: 

u  -  v  -  w  -  o  for  x  "  o 

w  *»  o  for  x  -  a  (43) 

w  -  o  for  y  »  o  and  y  -  b 
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The  total  potential  energy,  ir,  is  the  sum  of  the  strain  energy 
stored  in  the  skin,  Uw,  in  the  stringers,  Us,  in  the  frames,  Uf,  and  the 
potential  of  the  external  load,  Q,  and  is  written  as: 


»r  •  Uw  +  Us  +  Uf  +  0  (44) 

* 

The  strain  energy  in  the  skin  for  an  anisotropic  plate  with  A^g  - 

* 

A26  -  0  is  given  by 
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*  2  *  *  2  *  2 
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(45) 


*  2  *  *  2 
+  D22W,yy  +  4D26W,yyW,Xy  +  4D66w,xy)  dv 

where  ex,  ey  and  yxy  are  the  strain  components.  Commas  denote  differentiation 
with  respect  to  the  subscripted  variables. 


The  strains  are  expressed  in  terms  of  the  displacements  u,  v  and  w 
using  the  nonlinear  strain- displacement  relations: 
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In  the  derivations  that  follow,  the  coordinate  variables  x  and  y  are 
normalized  with  respect  to  their  respective  panel  dimensions.  The  normalized 
coordinates  (£,>?)  are  given  by 

£  -  *  ,  r,  -  X  (47) 
a  b 


The  strain  displacement  relations  can  then  be  rewritten  as: 
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In  Equation  (48)  and  hereafter,  the  subscripts  to  the  displacements  u,  v  and  w 
denote  differentiation  with  respect  to  the  subscript  variables. 


Substituting  the  strain- displacement  relations,  Equation  (48),  into 
Equation  (45)  the  strain  energy  stored  in  the  skin  becomes: 
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The  strain  energy  in  the  stringer  is 
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The  strain  energy  in  the  frame  is 


A  ~E 2  I_E~ 

Up - O  (1.  »7)  d»7  +  -£-£ 

Jo  2b3 


2b 


2 

(1 » »?)  d»? 


The  potential  of  the  external  loads  is 
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The  solution  method  employs  the  principle  of  minimum  potential 
energy.  In  applying  the  principle  of  minimum  potential  energy,  the 
displacement  components  are  assumed  to  be  functions  of  the  independent 
variables  4  and  V-  The  selected  functions  must  satisfy  the  displacement 
boundary  conditions  given  in  Equation  (43)  and  minimize  the  total  potential 
energy.  A  generalized  series  expression  for  the  displacement  functions  with 
unknown  coefficients  was  selected  for  the  present  analysis  and  these  are  as 


u  -  Anmfx  +  axa£ 

v  “  Bnmf2  +  blaS  (53) 

w  “  cnmf3  +  Dnmf4  +  w0f5> 

In  Equation  (53),  the  functions  fj  -  fi(x,y;n,m) ,  i  -  1,  2,  3,  4,  are 

arbitrarily  selected  admissible  functions.  The  expressions  such  as  Amfi  are 

shorthand  expression  for  a  double  series,  i.e., 

N  M 

Anmfi  *“  s  s  Anm  fl  (x,y;n,m) 
n-1  m=l 

The  coefficients  Am,  B^,  C^j,  D^,  a\  and  b^  are  unknown  coefficients  to  be 
determined  by  minimizing  the  total  potential  energy.  The  term  WQ  in  Equation 
(53)  is  the  initial  imperfection  at  the  panel  center.  The  function  f5"*f5(£,  ri) 
is  the  initial  imperfection  function  in  terms  of  the  lateral  displacement  and 
satisfies  the  displacement  boundary  conditions. 


Substituting  Equation  (53)  into  the  energy  expressions  and  using  the 
definitions  given  in  Appendix  A  for  the  individual  energy  integrals,  the 
energy  expressions  finally  become: 


U, 


w 


* 

A11 


2a 


11 

Anm  Apq  ^11  +  bal  Anm 


1  u  2 

Fl  +  &  al 


.  Ill  111  111 

+  ^nm  Cpq  Crs  Hi 3 3  +  2  A^  Cpq  Drs  H134  + 

Anm  Bpq  Drs  Hl44> 


bw  111  111  11 

+  __o  (Ajyjj  Cpq  H135  +  A^  Dpq  H145)  +  —  (al  Cpq  G33 
EL 


11  11  bw  11  11 

+  2al  Cjyjj  Dpq  G34  +  Dpq  G44)  +  — £  (al  C^  G35  +  a^  D^,  G45) 


bw2  111  bw2  11  ,  r  1111 

+  — ^  Anm  H155  +  — -  ai  G55  +  — —  Cnm  Cpq  Crs  Ctu  I3333 
2a2  2a  8a3  L 


1111  1111 
+  ^  ^nm  ^pq  Crs  Btu  -^3334  +  ^  ^nm  ^pq  Brs  Btu  *3344 
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1111  1111 
+  A  ^nm  Bpq  Drs  D^u  x3444  +  Bnm  Bpq  ^rs  ^tu  ^4444 

1111  1111 
+  4  wo  (cnm  cpq  crs  *3335  +  3  cnm  cpq  Drs  I3345 

1111  1111  1111 
+  3  Cj-jjjj  Dpq  Drs  I3445  +  Djuji  Dpq  Drs  I4445)  +  6  wQ  (C^,  Cpq  I3355 

1111  1111  1111 
+  2  Cm  Dpq  x3455  +  Bnm  Bpq  14455)  +  4  wo  (cnm  I3555 

1111  1111 
+  ^nm  x4555>  +  wo  I5555 

*  12  2 

+  a12  '  Anm  Bpq  G12  _  aal  Brun  F2 

. 


10  10  10 
+j|  (Anm  Gpq  G13  +  Anm  BpqG14  +  wo  Anm  C15) 


+  ^al  ^nm  ^3  +  al  Bnm  f4  +  wo  al  F5^ 

R 

122  122  122 
+  (Anm  cpq  crs  H133  +  2Anm  Cpq  Drs  H134  +  A^,,  Dpq  Drs  H344) 


w  122  122  22 

+  (Arun  Gpq  ^135  +  Anm  ^pq  ^145)  +  (al  C^,,  Cpq  G33 


22  22  aw  22 

+  2ax 

Grun  Bpq  G34  +  al  Bnm  Bpq  G44)  +  ~~  (al  Grun  G35 


22  w2  122  aw2  22 

+  al  ®run  g45)  +  ~~  Anm  ^155+  al  G55 
lb  lb 


,  211  211  211 
+  —  (Bran  Gpq  Grs  **233  +  2Bnm  Gpq  ^rs  **234  +  Bnm  **pq  ^rs  **244 


w  211  211  w2  211 

+  ^  (Bnm  Gpq  **235  +  Bnm  Dpq  **245 )  +  ^  Bnm  **255 


.  Oil  011  110 

+  - —  [Cnm  Cpq  Crs  H333  +  Cj^  Cpq  Drs  (2H334  +  H334)  (54  cont  d) 
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t 


A 


011  110  011 
+  cnm  Dpq  Drs  (H344  +  2H344)  +  Dnm  Dpq  Drs  H444 

011  110  011  101  110 
+  w0  Cnm  Cpq  (2H335  +  H335)  +  2w0  0^  Dpq  (H345  +  H345  +  345) 

011  110  011  101 
+  wQ  Dm,,  Dpq  (2H445  +  H445)  +  wQ  Cnm(H355  +  2H355) 

2  011  110  3  011 

+  w0  Djuj  (H455  +  2H455)  +  wQ  H555] 

,  r  1122  1122  1221 
+  —  Cnn,  Cpq  Crs  Ctu  I3333  +  2Cm  Cpq  Crs  Dtu  (I3334  +  I3334) 
4ab  L 


1122  1212  2211 
+  ^nm  ^pq  ^rs  ®tu  (13344  +  ^3344  +  13344) 

1122  2112  1122 
+  2CJUU  Dpq  Drs  Dtu  (I3444  +  I3444)  +  Dnm  Dpq  Drs  Dtu  J4444 

1122  1221  1122 
+  2w0  Cnn,  Cpq  Crs  (I3335  +  I3335)  +  2w0  Cpq  Drs  (I3345 

1212  1221  2211  ^  ,1122  2112 
+  213345  +  213345  +  13345)  +  2w0  Cryn  Dpq  Drs  (213445  +  !3445 


1221  2211  ,  U22  1221 

+  13445  +  2I3445)  +  2w0  Djjjji  Dpq  Drs  (I4445  +  I4445) 


1122  1212  2211  1122 
+  wo  cnm  cpq  (13355  +  4I3355  +  !3355)  +  2w0  Cnm  Dpq  (I3455 


1212  2112  2211  2  1122  1212  2211 
+  2i3455  +  213455  +  13455)  +  wo  ^nm  ^pq  (14455  +  ^4455  +  14455 

3  1122  2112  3  1122  2112  4  1122  ' 

+  2wc  (I3555  +  13555)  +  2wo  Dnm  (14555  +  x4555)  +  wo  !5555  . 


* 

+  A22 


a 

2b 


22 

®nm  ®pq  G22 


20 


20 


+  “  (®nm 
K 


Cpq  O23  +  Bnjji  Dpq  G24 


20  00  00  00 
+  w0  Bnm  G25)  +  (Cnm  Cpq  G33  +  2Cnm  Dpq  G34  +  Dnm  Dpq  G44> 
2R2 


abw  00  00  abw2  00  222 

+  (Cl™  G35  +  D“"  G45>  +  ^  °55  +  P<1  IS  ”54  Cont'd) 
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222  222  aw  222 

+  2Bnm  Cpq  Drs  H234  +  Bnm  Dpq  Drs  H244)  +  (Bnm  Gpq  H235 

222  aw2  222  r  „  ,022 

+  Bnm  Dpq  H245>  +  ^  Bma  Hz55  +  ^  °Pq  ^  H333 


022  220  022  202 
+  cnm  cpq  Drs  (2H334  +  H334)  +  cnm  Dpq  Drs  (H344  +  2H344) 


022 

+  Onm  Opq  Drs  H444 


1  aw 
+  — fi  L 

J  2bR 


cnm  cpq  (2H335  +  2H335  +  H335> 


022  202  220  022  220 
+  2Cnm  Dpq  (H345  +  H345  +  H345)  +  Dnm  Dpq  (2H445  +  H445) 


aw2  [  022  202  022  202  1  aw3  022 

+  —  L  Cnm  CH355  +  2H355)  +  Dnm  (H455  +  2H455>  J  +  ^  H555 


2bR 


+  _a_  [ 


2222 


2222 


8b3 


cnm  cpq  crs  ctu  I3333  +  ^cnm  cpq  crs  Dtu  13334 


2222 


2222 


+  6Cnm  cpq  Drs  Dtu  *3344  +  ^cnm  Dpq  Drs  Dtu  I3444 


2222 


2222 


+  Dnm  Dpq  Drs  Dtu  *4444  +  Aw0  (cnm  cpq  crs  I3335 

2222  2222  2222 
+  3cnm  Cpq  Drs  I3345  +  3Cnm  Dpq  Drs  I3445  +  Dnm  Dpq  Drs  *4445 

2  2222  2222  2222 
+  6w0  (Cnm  cpq  13355  +  2cnm  Dpq  I3455  +  Dnm  Dpq  I4455) 

3  2222  2222  2222 
+  4w0  (Cnm  13555  +  Dnm  *4555)  +  wo  I5555  J 


* 

+  a66 


22  11  1  .  2 
r*  Anm  Apq  Gn  +  ^  Bnm  »pq  G22  +  bbl  Bnm  F2  +  ^  bl 
2b  z 


21  2  1  r  212 

Anm  Bpq  G12  +  abl  Anm  F1  +  £  I  Anm  cpq  crs  H133 
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+ 
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+ 


212  221  212 
Anm  cpq  Drs  (H134  +  H134)  +  Anm  Dpq  Drs  H144 

221  212  212  221 
wo  Anm  cpq  (**135  +  H135  +  wo  Anm  Dpq  <H145  +  H145) 

1  ,  [  112  112  121 

wo  Arun  **155  J  +  —  *-  Bm  Cpq  Crs  h233  +  Bnm  cpq  Drs  (H234  +  **234> 
& 

112  121  112 
Bnm  Bpr  '*s  **244  +  wo  Bnm  Cpq  (H235  +  **235> 

112  121  2  112  I  12 

wo  Bnm  Dpq  (**245  +  H245)  +  wo  Bnm  H255  J  +  bl  cnm  cpq  633 

12  12 
bl  Cnm  Bpq  (C34)  bl  Bnm  Bpq  ®44 

r  21  12  21  12  I  2  12 

wQ  bl  Cnm  (G35)  +  G35)  +  bx  0^(645  +  G45)J  +  w0  bX  G55 

r  1122  1122  1221 
~r  C^n  Cpq  Crs  Ctu  I3333  +  2Cnm  Cpq  Crs  Ddtu  II3334  +  13334) 
2ab  L 

1122  1212  2211 

Cnm  Cpq  Drs  D^u  (x3344  +  ^x3344  +  x3344) 

1122  2112  1122 
2Cnm  Bpq  Brs  Btu  (*3444  +  x3444)  +  Bnm  Bpq  Brs  Btu  x4444 

1122  1221 

2w0  Cjjjj  Cpq  Crs  (I3335  +  13335) 

1122  1212  1221  2211 
2w0  Cm,,  Cpq  Drs  (I3345  +  2I3345  +  2I3345  +  I3345) 

1122  2112  1221  2211 
2w0  Cnm  Ddpq  Drs  (2I3445  +  13445  +  I3445  +  2I3445 

1122  1221 


2w0  Bnm  Bpq  Brs  (x4445  +  x4445) 

2  1122  1212  2211 

wo  cnm  Cpq  (13355  +  4I3355  +  I3355> 

2  1122  1212  2112  2211 
2wo  cnm  Dpq  (I3455  +  2I3455  +  2I3455  +  x3455) 

2  1122  1212  2211  3 

wQ  Dnm  Dpq  <x4455  +  4I4455  +  I4455>  +  2wo  cnm 
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3  1122  2112  4  1122 

+  2w0  Dm  (I4555  +  I4555  +  w0  I5555 


bD*  33  33  33  33 

+  — il  (Cju,,  Cpq  G33  +  20^  Dpq  G34  +  Dm,,  Dpq  G44  +  2w0  G35 

2a3 


Jpq 

33  2  33 


pq 

34 


+  2w0  Djyj  G45  +  wQ  G55)  +  —j£  Cnm  Cpq  G33 

34  43  34  43  34 

+  Gnm  ^pq  (G34  +  G34)  +  Dnm  ^pq  G44  +  wo  Gnm  (G35  +  G35) 


+  wo  ®nm  (G45 


34  43  2  34  1  2D*  I 

45  +  G45>  +  w0  G55  J  +  —la  L 


35 

Gnm  Gpq  G33 


35  53  35  35  53 

+  Gnm  ®pq  (G34  +  G34)  +  ®nm  ^pq  G44  +  wo  Gnm  (G35  +  G35) 


35 

+  wo  ®nm  (G45 


53  2  35  1  aD*  f 

+  G45)  +  w0  G55  I  +  —M  L 

2b3 


44 

Gnm  Gpq  G33 


44  44  44  44 

+  2Cnm  Dpq  G34  +  Djyn  Dpq  G44  +  2w0  G35  +  2w0  D^  G45 


44 

+  wo  G55 


]+2Ss[ 

J  b2 


45  45  54 

Gnm  Gpq  G33  +  Gnm  Dpq  (G34  +  G34) 


45 


45  54 


45  54 


■*"  Dnm  Dpq  G44  +  wQ  (G35  +  G35)  +  w0  Dj-m,  (G45  +  G45) 


2  45  ]  2D*  [ 

+  wo  g55  J  +  L 


55  55 

Gnm  Gpq  G33  +  20,^  Dpq  G34 


55  55  55  55 

+  Djjjj  D(jpq  G44  +  2w0  G35  +  2w0  Dj^  G45  +  w0  G55 


Ue 


A  E  [  11  12  2 

L  ^nm  ^pq  ^11  +  2a  a]_  +  a^a 


s  s 
2a 


I  E 
s  s 

2a3 


33 

^nm  ®pq  ^22 


A..E.  22  I-E-  44 

“  nc  ®nra  ®pq  ^22  +  3  ^pq 


2b 
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n  -  -  Nxx  [  Am  JX  +  aia  ]  b 
-  Nxy  [  Bnm  J2  +  bla  ]  b 


The  total  potential  energy  is  minimized  with  respect  to  the  unknown 
coefficients.  The  minimization  process  yields  a  system  of  nonlinear  algebraic 
equations.  Details  of  these  algebraic  equations  are  given  in  Appendix  A. 
These  equations  can  be  expressed  in  the  following  form: 


where  the  subscript  L  denotes  the  linear  terms  of  the  partial  derivative  of 
the  total  potential  energy  with  respect  to  a  particular  unknown  coefficient 
(A) ,  subscript  N  denotes  the  nonlinear  terms  and  C  represents  the  terms  that 
are  independent  of  the  unknown  coefficients. 


An  alternate  approach  to  minimize  the  total  potential  energy  is  to 
directly  evaluate  the  energy  expressions  given  by  Equations  (54)  through  (57). 
In  this  approach,  the  values  of  the  unknown  coefficients  are  initially 
assumed.  The  value  of  the  total  potential  energy  is  then  directly  evaluated 
using  the  assumed  values  of  the  coefficients.  The  values  of  the  coefficients 
are  systematically  varied  and  the  value  of  the  corresponding  total  potential 
energy  is  evaluated  iteratively  until  the  minimum  potential  energy  is 
approximately  achieved.  The  values  of  the  unknown  coefficients  are  then 
substituted  into  Equation  (58)  to  verify  the  numerical  accuracy.  This 
procedure  is  adopted  as  an  alternate  approach  because  of  a  convergence  problem 
encountered  in  solving  the  nonlinear  system,  Equation  58,  directly. 

The  numerical  solution  procedures  are  discussed  in  Sections  3.2  and 
3.3.  Section  3.2  details  a  single  mode  solution,  in  which  only  one  buckling 
mode  is  selected  for  analysis,  i.e.,  only  one  term  for  each  of  the 
displacement  functions  fj_,  i  -  1,  2,  3,  4,  is  used.  In  this  case,  the  system 
of  Equations  58  is  reduced  to  six  equations  and  direct  solution  of  Equation  58 
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presents  little  difficulty.  The  multi-mode  solution  is  discussed  in  Section 
3.3.  Because  of  the  interaction  between  buckling  modes,  numerical  solution  of 
Equation  58  is  not  always  possible.  A  combination  of  the  two  approaches 
discussed  earlier  is  used  to  determine  the  unknown  coefficients. 

In  the  actual  numerical  solution,  the  displacement  functions  hav'-  to 
be  specifically  defined.  The  functions  that  satisfy  the  boundary  conditions 
and  describe  the  displacement  behavior  observed  in  postbuckled  panel  tests  are 
as  follows: 

f^(£,  «| ;  n,  m)  «  sin  II£  cos  (nll£  -  mlty) 
f2(£.  ri;  n,  m)  -  (1  -  cos  II£)  cos  (nli£  -  mllr?) 

f3<£,  *i‘>  n,  m)  -  sin  nH£  sin  mllrr  (59) 

f4<£,  i||  n,  a)  ■  sin  II£  sin  Tlri  sin  (nll£  -  mllr;) 

f5<f ,  t\)  -  sin  n£  sin  Ify 

The  individual  energy  integrals  written  in  terms  of  the  above 
functions  are  evaluated  in  closed  form.  Detailed  expressions  of  all  the 
energy  integrals  are  given  in  Appendix  A. 

3.2  SINGLE  MODE  ANALYSIS 

The  number  of  nonlinear  algebraic  equations  in  the  system  given  by 
Equations  (58)  depend  on  the  number  of  buckling  modes  used  in  the  analysis. 
The  number  of  equations  can  be  calculated  from  the  relation  4NM+2,  where  N  is 
the  number  of  buckling  modes  in  the  x-direction  and  M  is  the  number  of 
buckling  modes  in  the  y-direction.  As  N  and  M  increase,  the  number  of 

nonlinear  terms  on  the  right  hand  side  of  Equation  58  also  significantly 

increases  (see  Equations  A-9  through  A-14  in  Appendix  A).  A  large  number  of 

nonlinear  terms  present  numerical  difficulty  in  solving  Equation  (58).  On  the 
other  hand,  although  the  postbuckling  behavior  of  a  stiffened  panel  is  mixed¬ 
mode  behavior  in  general,  the  displacement  response  is  dominated  by  a  single 
buckling  mode.  Therefore,  if  the  dominant  buckling  mode  is  known,  the 
postbuckling  behavior  can  be  accurately  described  using  a  single-mode 
analysis. 
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Let  n  be  the  selected  buckling  mode  in  the  x- direction  and  m  the 
mode  in  the  y-direction.  Then  the  number  of  equations  in  the  nonlinear  system, 
Equation  (.58)  is  reduced  to  six.  Since  there  is  no  mode  interaction,  the 
energy  integrals  G,  H,  I,  K  in  Equations  A- 9  through  A- 14  become 


af)  afi 

Gij ,  nmpq  "  Gij ,  nmnm 


afo  aM 

Mijk,nmpqrs  “  ni j k , nmnmnm 


af}-yS  af3-yS 

*ijkl,  nmpqrstu  *  *ijkl,  nmnmnmnm 


^aa.nmpq  *  *aa*  nmnm 


(60) 


In  addition,  for  the  case  of  no  initial  imperfection,  i.e.,  wQ  -  o,  the  total 
number  of  integrals  involved  in  Equations  58  is  reduced  to  108.  The  reduced 
system  of  nonlinear  equations  can  be  solved  with  "  :y  high  accuracy  by  an 
iterative  technique  using  the  method  of  successive  linearization.  In  this 
method,  each  of  the  unknown  coefficient,  A^,  B^,  C^,  D^,  a^  and  b^  is 
assigned  an  initial  value  and  substituted  into  the  right -hand- side  of  Equation 
58.  Equation  58  now  becomes  a  system  of  linear  algebraic  equations  and  can  be 
easily  solved  for  the  new  values  of  the  unknown  coefficients.  Using  the  new 
set  of  coefficients  as  initial  values,  another  set  of  improved  coefficients 
can  be  obtained  by  solving  the  linearized  system.  This  procedure  is  continued 
until  the  solution  converges  within  a  desired  limit.  In  the  actual  solution, 
only  the  initial  values  of  the  coefficients  at  the  first  load  level  need  to  be 
assigned.  At  higher  load  levels,  the  initial  values  are  obtained  by 
extrapolating  the  converged  solutions  of  the  preceding  load  levels  to  reduce 
the  number  of  iterations. 


The  results  of  the  single -mode  solution  for  an  example  problem  are 
presented  below.  The  panel  geometry  and  the  material  properties  used  in  this 
problem  are: 

Panel  Length  a  -  17.5  in. 

Panel  Width  b  -  10.0  in. 

Panel  Radius  R  -  45.0  in. 
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★ 

An  -  562  kips/in 
* 

A22  **  582  kips/in 
* 

DU  -  121  lb/in 
* 

■ie 

A12  -  174.7  kips/in 
* 

Agg  -  225  kips/in 
* 

D12  -  36.4  lb/in 
* 

* 

°16  “  d26  -  0 

D22  “  121  lb/in 

Dgg  -  46.87  lb/ in 

As  -  0.11  in2 

Es  -  10  x  10^  psi 

Is  -  0.00615  in4 

Af  -  0.5  in2 

Ef  -  10  x  106  psi 

If  -  5.0  in4 

In  this 

example,  the  frame  is  assumed 

to  be  very  rigid,  this 

simulated  by  using  a  relatively  large  moment  of  inertia  (5.0  in4)  as  compared 
to  that  of  the  stringer  (0.00615  in4).  The  ratio  of  the  compression  to  shear 
load  (Nx/NXy)  used  in  the  example  is  -1.0. 

The  results  of  the  single  mode  analysis  for  this  example  are 
illustrated  in  Figures  8  through  15.  These  results  include  solutions  for  m  =  1 
with  n  rarging  from  1  to  6.  The  value  of  a^  as  a  function  of  the  total 
compression  load  (Pxx  -  bNx)  is  shown  in  Figure  8.  The  parameter  a\  is  an  in¬ 
plane  displacement  parameter.  The  value  of  af  is  essentially  the  axial  strain 
in  the  x-direction  due  to  end- shortening  and  is  obtained  from  the  first  of 
Equations  53.  As  shown  in  Figure  8,  a^  varies  linearly  with  applied  load  for 
n  -  1.  For  n  t  2,  the  end- shortening  parameter  becomes  approximately  bilinear 
with  applied  load.  The  load  level  where  the  slope  of  the  end- shortening  curve 
changes  signifies  initial  buckling  of  the  skin.  The  figure  shows  that  the 
initial  load  is  lower  for  some  of  the  higher  buckling  modes  (larger  value  of 
n) .  The  initial  buckling  load  for  n  -  2  is  approximately  5000  lbs .  and 
reduces  to  2800  lbs  for  n  -  6.  These  results  indicate  that  with  m  -  1 ,  the 
first  two  buckling  modes  (n  -  1  and  n  -  2)  are  not  likely  to  dominate  the 
postbuckling  behavior  of  the  panel.  This  conclusion  is  more  readily  evident 
from  the  results  of  other  parameters  discussed  below. 

Figure  9  shows  the  relationship  between  the  shear-displacement 
parameter  bf  and  the  total  applied  compression.  The  parameter  bf  is  the 
dominant  term  in  the  y-direction  displacement  ( shear- displacement) .  The  value 
of  b^a  is  the  average  shear -displacement  at  the  loading  end  of  the  panel  (x  - 
a) .  The  results  shown  in  Figure  9  indicate  that  the  value  of  bf 
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Figure  9.  Shear  End-Displacement  Parameter  as  a  Function  of  the 
Applied  Compression  Load. 
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IN-PLANE  DISPLACEMENT  PARAMETER,  Anmx10 


Figure  10.  Displacement  Coefficient  Anm  as  a  Function  of  the  Applied 
Compression  Load. 
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OUT-OF-PLANE  DISPLACEMENT  PARAMETER, 


Figure  12.  Displacement  Coefficient  C^  as  a  Function  of  the  Applied 
Compression  Load. 
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OUT-OF-PLANE  DISPLACEMENT  PARAMETER 


TOTAL  APPLIED  COMPRESSION,  -Pxx/103  (lb) 

Figure  13.  Displacement  Coefficient  Dnra  as  a  Function  of  the  Applied 
Compression  Load, 
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Figure  14.  Maximum,  Minimum  and  Panel  Center  Out- of- Plane  Displacements 
as  Functions  of  the  Applied  Compression  Load. 
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Figure  15.  Out-of-Plane  Displacement  Contours  for  the  Buckling  Mode 
n  *=  3,  m  =  1,  at  an  Applied  Compression  Load  of  8500  lbs. 


varies  linearly  with  the  applied  compression  for  n  -  1.  For  n  2:  2,  the  bi 
curves  are  slightly  nonlinear.  The  values  of  are  approximately  equal  for 
n  -  3  to  6. 

The  in-plane  displacement  parameters  A^,  and  b^  are  shown  in 
Figures  10  and  11.  The  coefficient  A^,  shown  in  Figure  10  is  a  measure  of  the 
in-plane/out-of-plane  displacement  interaction  in  the  x-direction.  As  shown 
in  the  figure,  for  n  £  3,  this  parameter  remains  approximately  zero  when  the 
applied  load  is  below  the  initial  buckling  load.  At  higher  loads  the  absolute 
value  of  the  parameter  increases  with  applied  load.  For  n  -  1  and  2,  the 
absolute  value  of  A^  continuously  increases  with  applied  load.  These  results 
again  indicate  that  the  panel  analyzed  is  not  likely  to  deform  into  the  first 
two  modes  (n  -  1,  2).  Figure  11  shows  a  similar  behavior  for  the  coefficient 
Bnm>  the  in-plane/out-of-plane  displacement  interaction  coefficient  for  the  v- 
displacement. 

Figures  12  and  13  show  the  variation  of  the  out-of-plane 
displacement  coefficients  and  Dm.  The  trends  for  these  displacement 
parameters  are  similar  to  that  of  A^  and  B^.  However,  they  are  one  order  of 
magnitude  higher  than  Anm  and  Bnra. 

Figure  14  shows  the  maximum,  minimum  and  panel-center  out-of-plane 
displacement  variations  with  the  total  applied  compression  load  for  the 
buckling  mode  n  -  1,  m  -  3.  The  panel  center  is  at  (8.75,  5).  As  shown  in 
the  figure,  the  out-of -plane  displacement  at  the  panel-center  remains 
relatively  small  as  the  load  increases.  This  is  because  the  center  is  in  the 
vicinity  of  a  nodal  line  for  the  assumed  buckling  mode.  The  maximum  (outward) 
displacement  occurs  at  location  (6.25,  5).  Figure  14  shows  that  the 
displacement  at  this  point  remains  at  approximately  zero  below  a  compression 
load  of  3500  lbs.  Above  3500  lbs.,  the  displacement  increases  rapidly  as  the 
load  approaches  4000  lbs.  Beyond  4000  lbs.,  the  displacement  increases  with 
load  at  a  relatively  slower  rate.  The  minimum  (maximum  inward)  displacement 
occurs  at  the  location  (11.25,  5).  The  out-of-plane  displacement  at  this 
point  varies  with  applied  compression  load  in  a  manner  similar  to  that  of  the 
maximum  displacement  but  in  the  opposite  direction. 
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The  out-of-plane  displacement  contours  for  the  buckling  mode  n  -  3, 
u  «  1  at  applied  compression  load  of  8500  lbs.  are  shown  in  Figure  15.  The 
figure  shows  that  there  are  two  major  buckles  in  the  center  portion  of  the 
panel.  These  buckles  are  oriented  at  an  angle  of  approximately  60°  from  the 
x-direction.  The  buckle  near  the  fixed  end  (left  buckle)  deforms  outward 
(positive  displacement) .  The  buckle  near  the  loading  end  (right  buckle) 
deforms  inward  (negative  displacement).  In  addition  to  these  major  buckles, 
two  minor  adjacent  buckles  also  develop.  These  are  shown  at  the  left-upper 
and  right- lower  corner  in  Figure  15. 

The  results  shown  above  indicate  that  the  single-mode  analysis 
developed  here  can  be  used  to  investigate  the  postbuckling  behavior  of  a  panel 
in  several  ways.  First,  a  parametric  study  can  be  carried  out  to  determine 
the  possible  buckling  modes  of  a  panel.  Once  a  dominating  mode  is  selected, 
the  single-mode  analysis  can  be  used  to  describe  the  approximate  postbuckling 
displacement  pattern  of  the  panel.  This  analysis  provides  the  approximate 
displacement  and  stress  (strain)  fields  of  the  panel  at  different  load  levels. 
This  information  can  then  be  used  for  a  fatigue  or  failure  analysis  of  the 
panel . 

3.3  MULTI -MODE  ANALYSIS 

In  the  single  mode  analysis,  the  solution  is  obtained  for  a 
preselected  buckling  mode.  The  implicit  assumption  in  the  single  mode  analysis 
is  that  the  panel  deforms  in  a  fixed  mode  as  the  applied  load  increases  in  the 
postbuckling  regime.  This  assumption  has  been  adopted  in  a  number  of  analyses 
in  the  literature.  However,  experimental  data  Indicated  that  buckling  mode 
shape  may  change  as  the  applied  load  increases.  In  order  to  simulate  the 
change  of  buckling  mode  in  the  postbuckling  regime,  a  multi-mode  analysis 
method  is  needed. 

The  multi-mode  analysis  method  developed  is  based  on  the  strain 
energy  method  discussed  in  3.1.  However,  a  general  solution  of  Equation  58  is 
prohibitive  because  of  the  excessive  computing  resource  requirements  and 
convergence  problems  in  the  numerical  solution.  The  number  of  nonlinear 
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equations  in  Equation  58  is  4NM+2  for  NxM  inodes  used  in  the  general  analysis. 
The  number  of  energy  integrals  increases  rapidly  with  N  and  M.  Table  2  shows 
the  total  number  of  energy  integrals  required  for  different  values  of  N  and  M. 
The  numbers  shown  in  Table  2  include  the  initial  Imperfection  terms.  These 
numbers  suggest  that  selected  buckling  modes  up  to  4  x  1  would  be  more 
practical. 

In  addition  to  the  number  of  energy  integrals,  the  number  of 
equations  also  limits  the  number  of  buckling  modes  used  in  the  analysis.  This 
is  because  the  interaction  of  buckling  modes  results  in  convergence  problems 
in  solving  the  nonlinear  equations.  In  the  multi-mode  solution,  the  two 
numerical  procedures  discussed  in  3.1  for  either  solving  Equations  58  or 
directly  minimizing  the  total  potential  energy  are  used  alternately  to  avoid 
a  convergence  problem. 

The  analysis  procedure  for  the  multi-mode  solution  was  coded  in  a 
Fortran  computer  program  PACL  (Postbuckling  Analysis  for  Combined  Loads) . 
Details  of  the  computer  program  are  documented  in  a  separate  volume  of  this 
report  -  Automated  Data  Systems  Documentation  (Reference  14) . 
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Table  2.  Ntiraber  of  Integrals  Required. 


NUMBER  OF 

N  TERMS 

NUMBER  OF 

M  TERMS 

NUMBER  OF 
CONSTANTS 

1 

1 

230 

2 

1 

1222 

3 

1 

3976 

4 

1 

9944 

4 

2 

110,956 

5 

1 

21,034 

6 

1 

39,610 

6 

2 

498,112 

7 

1 

68,492 

8 

1 

110,956 

8 

4 

21,725,092 
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SECTION  4 

ANALYSIS  AND  TEST  RESULTS  CORRELATION 


4.1  Introduction 

The  static  and  fatigue  data  presented  in  Reference  15  were 
*  analyzed  to  correlate  the  measured  initial  buckling  load,  ultimate  strength 

and  failure  mode  with  predictions  from  the  semi -empirical  analysis 
r  methodology.  The  measured  strains,  buckling  mode  shape  changes  and  panel 

stiffness  changes  were  compared  with  predictions  from  the  energy  method  based 
analysis  presented  in  Section  3.  The  fatigue  life  data  were  utilized  to 
establish  S-N  curves  for  metal  and  composite  panels.  These  results  are 
discussed  in  the  following  paragraphs. 

4.2  Initial  Buckling  Under  Combined  Loads 

Initial  buckling  predictions  for  the  metal  panels  were  based  on 
the  semi-empirical  expressions  (Reference  6)  given  in  Equations  1  and  3.  A 
comparison  of  the  buckling  load  data  obtained  in  this  program  for  pure  shear 
or  pure  compression  loading,  with  the  predictions  is  shown  in  Figure  16.  As 
seen  in  the  figure,  the  semi -empirical  predictions  are  conservative  by  as  much 
as  30  percent.  This  observation  is  consistent  with  the  metal  panel  data 
generated  in  Reference  1. 

Under  combined  loading,  the  buckling  loads  were  predicted  using  a 
parabolic  interaction.  A  comparison  of  the  predictions  and  the  test  data  is 
shown  in  Figure  17.  Since  the  pure  shear  and  pure  compression  buckling  loads 
were  higher  than  predicted,  the  buckling  loads  under  combined  compression  and 
shear  are  underestimated  by  Equation  4.  A  true  comparison,  however,  can  be 
seen  in  Figure  18  where  the  interaction  is  shown  in  terms  of  the  buckling  load 
4  ratios  Rc  and  Rs.  The  comparison  in  Figure  18  shows  excellent  correlation 

between  the  predictions  and  the  test  data  provided  the  scatter  in  the  measured 
values  for  the  different  panels  is  accounted  for  by  normalization.  Reference 
16  data,  used  to  establish  the  parabolic  interaction  equation, are  also  shown 
to  illustrate  the  consistency  of  the  data  obtained  in  the  present  program. 
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COMPRESSION  BUCKUNG  LOAD  N,  ,  Ibs/in 


Figure  18  shows  that  the  buckling  interaction  equation  used  provides  a  lower 
bound  and  is  thus  somewhat  conservative  for  buckling  load  predictions  for 
metal  panels  under  combined  shear  and  compression  loading. 

The  correlation  between  initial  buckling  load  predictions  and  test 
data  for  composite  panels  is  shown  in  Figure  19.  In  the  case  of  composite 
panels,  the  pure  shear,  pure  compression  and  the  combined  loading  initial 
buckling  predictions  were  based  on  program  SS8  (Reference  8).  In  Figure  19, 
the  parabolic  and  linear  interaction  curves  are  also  shown  for  comparison. 
The  linear  interaction  expression  provides  a  lower  bound  for  the  test  data. 
For  preliminary  design  purposes,  use  of  the  linear  interaction  is  more 
appropriate  for  composite  panels.  Figure  20  shows  that  the  test  data  are 
bounded  by  the  linear  and  a  fourth  power  (i.e.,  a  -  1  and  a  -  4  in  the  expres- 
a 

sion  Rc  +  Rs  -  1) . interaction  rule  (Reference  17). 


4.3  ULTIMATE  STRENGTH  UNDER  COMBINED  LOADS 


The  ultimate  strength  of  metai  and  composite  panels  was  predicted 
using  the  methodology  given  in  Section  2.  The  strength  predictions  were 
plotted  as  failure  envelopes  and  are  shown  in  Figures  21  and  22  for  metal  and 
composite  panels,  respectively.  The  only  change  in  the  strength  prediction 
methodology  made  after  comparison  with  test  data  was  in  the  stiffener 
crippling  interaction  equation  under  combined  loading.  Originally  (Equation 
41)  the  following  criterion  was  adopted  for  stiffener  crippling: 
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However,  the  test  data  for  both  metal  and  composite  panels  show  better 
correlation  with  a  linear  interaction,  i.e., 
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As  shown  in  Figure  21  the  metal  panel  strength  data  show  excellent  correlation 
with  predictions.  The  failure  mode,  however,  was  permanent  set  in  the  skin. 
The  close  agreement  of  the  test  data  with  stiffener  crippling  predictions 
leads  to  the  conclusion  that  the  skin  creasing  was  precipitated  by  stiffener 
crippling. 


The  composit  s  panel  test  data  also  show  good  agreement  with  the 
linear  interaction  stiffener  crippling  prediction.  There  are  two  exceptions, 
however,  in  panels  GR-1  and  GR-2.  The  low  failure  loads  obtained  for  these 
panels  are  plausible  since  these  two  early  panels  showed  some  load 
introduction  problems  during  the  static  tests.  Specifically,  the  panel  load 
introduction  area  skin  thickness  was  the  same  as  the  test  section  skin 

thickness.  Due  to  load  introduction  eccentricities,  the  skin  in  the  load 
introduction  area  buckled  before  the  panel  ultimate  load  was  reached.  Thus, 
the  two  panels  were  not  subjected  to  a  uniform  axial  compression  load  and, 
therefore,  showed  failure  loads  slightly  lower  than  the  predictions.  In  all 
other  panels  the  load  introduction  region  thickness  was  increased  by 

secondarily  bonding  fiberglass  laminates.  Thus,  the  semi -empirical  design 
method  as  given  in  Section  2  with  Equation  61  replacing  Equation  41  can  be 
used  for  designing  curved  composite  panels  under  uniaxial  compression  and 
shear  loads. 

4.4  FATIGUE  LIFE  UNDER  COMBINED  LOADS 

The  fatigue  test  data  for  metal  and  composite  panels  are  fully 
documented  in  Reference  15.  A  summary  of  the  fatigue  failure  modes  for  the 
metal  panels  under  compression  dominated  constant  amplitude  loading  i.e. 
(Nx)max/(Nxy)max  "  2  with  R-10  for  compression  and  R—  1  for  shear*,  is  shown 

in  Table  3.  For  shear  dominated  constant  amplitude  fatigue  loading  i.e. 

(Nx)max/(Nxy)max  “0.5  with  R-10  for  compression  and  R--1  for  shear,  the  test 
results  and  failure  modes  are  shown  in  Table  4.  The  basic  fatigue  failure 
mode  in  the  metal  panels  under  compression  dominated  loading  was  crack 
initiation  in  the  skin  adjacent  to  the  stiffener  flange  and  subsequent 
propagation  along  the  loading  direction.  The  crack  initially  propagated  along 


♦Note  that  due  to  differences  in  R-ratios,  (Nx)max  and  (Ny)aax  do  not  occur 
simultaneously . 
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Table  3 .  Fatigue  Failure  Modes  of  Aluminum  Panels  Under 
Compression  Dominated  Loading. 
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the  stiffener  direction.  After  a  certain  length,  the  crack  branched  and  grew 
toward  the  centerline  of  the  bay  in  the  diagonal  direction.  The  crack 
initiation  life  of  the  panels  was  approximately  60  percent  of  the  total  number 
of  cycles  required  to  tear  one  skin  bay. 

Under  shear  dominated  loading,  cracks  initiated  at  the  edges  of 
fastener  holes  in  the  skin.  The  subsequent  crack  growth  pattern  and  the  crack 
initiation  life  relative  to  the  total  number  of  cycles  required  to  tear  one 
bay  were  similar  to  those  obtained  in  the  compression  dominated  tests. 

In  practical  aircraft  structures  curved  panels  are  most  commonly 
used  in  pressurized  fuselage  structures.  Therefore,  the  crack  initiation  life 
was  defined  as  the  fatigue  life  of  the  metal  panels.  It  should  be  noted, 
however,  from  the  test  data  that  postbuckled  metal  panels  retain  a  significant 
percentage  of  their  static  strength  even  after  the  loss  of  skin  due  to 
cracking  under  diagonal  tension  stresses. 

The  measured  crack  initiation  life  of  the  metal  panels  was  used  to 
generate  the  S-N  curve  shown  in  Figure  23.  In  this  figure,  the  typical 

fatigue  failure  mode  is  also  illustrated.  The  limited  data  generated  in  this 
program  indicate  steep  S-N  curves  for  metal  panels  operating  in  the 
postbuckling  range.  Secondly,  the  metal  panels  can  sustain  static  loads  of 
approximately  2.5  times  the  average  buckling  load  but  in  actual  structures 
their  capability  would  be  limited  to  1.25  times  the  buckling  loads  due  to 
fatigue  considerations. 

Fatigue  test  results  for  the  composite  panels  are  summarized  in 

Table  5.  The  two  panels  tested  at  (Nx)niax/(NXy)max  -  2  experienced  no  fatigue 

failure  after  100,000  cycles  of  constant  amplitude  loading.  Residual  static 

«  strength  tests  on  these  panels  indicated  no  strength  reduction  (See  Figure 

22).  The  static  failure  mode  was  primarily  skin/stiffener  separation.  Panels 
*•  under  constant  amplitude  shear  dominated  loads  i.e.  (Nx)raax/(NXy)raax  -  0.5, 

failed  under  fatigue  cycling.  The  fatigue  failure  mode  in  these  panels  (GR-7 
and  GR-8)  was  skin  stiffener  separation  at  stiffener  and  ring  intersection 
accompanied  by  local  skin  rupture,  Thus,  the  composite  panels  appear  to  be 
more  sensitive  in  fatigue  to  shear  dominated  loading. 
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SKIN  RUPTURE - STIFFENER/SKIN  DISBONDS 
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NUMBER  OF  CYCLES 


Figure  24  shows  a  plot  of  the  number  of  fatigue  cycles  sustained  by 
the  composite  panels  versus  the  applied  loads.  From  four  data  poircs  in 
Figure  24  a  fatigue  threshold  was  estimated  to  be  approximately  80  percent  of 
the  static  strength.  The  fatigue  advantage  of  composite  panels  relative  to 
metal  panels  is  readily  apparent  from  Figure  24  in  that  the  composite  panels 
could  be  utilized  up  to  200  percent  of  their  initial  buckling  load  for  shear 
dominated  loading  as  opposed  to  125  percent  for  metal  panels.  The 
postbuokling  range  for  composite  panels  under  compression  dominated  loading 
could  be  possibly  higher. 
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NUMBER  OF  CYCLES 


SECTION  5 


CONCLUSIONS 


The  significant  conclusions  from  this  program  are  summarized  in  the 
following  paragraphs. 

5.1  SEMI -EMPIRICAL  DESIGN  METHODOLOGY  FOR  POSTBUCKLED  PANELS  UNDER 

COMBINED  LOADS 

1.  The  semi-empirical  static  design  methodology  developed  in 
Reference  1  for  postbuckled  composite  and  metal  panels  under 
pure  shear  or  pure  compression  loading  was  extended  to  panels 
under  combined  uniaxial  compression  and  shear  loads. 

2.  The  methodology  was  coded  in  a  computer  program  (PBUKL)  for 
rapid  iterative  design  of  composite  and  metal  panels. 

3.  Experimental  verification  data  were  used  to  develop  a  new 

criteria  to  predict  the  effect  of  shear  and  compression  load 
interaction  on  composite  panel  skin  buckling.  A  linear 

interaction,  although  conservative,  seems  more  appropriate  for 
the  design  of  composite  panels  as  opposed  to  the  well 
established  parabolic  interaction  rule  for  metal  panels .: 

4.  The  test  data  showed  that  for  both  composite  and  metal  panels  a 
linear  interaction  rule  for  stiffener  crippling  prediction 
yields  better  correlation  than  a  non-linear  interaction  rule. 

5.  Ultimate  panel  strength  predictions  based  on  the  semi-empirical 
analysis  for  composite  and  metal  panels  were  found  to  be  very 
accurate  and  well  suited  for  design  purposes. 

6.  Stiffener  and  skin  separation  in  composite  panels  was  the 
observed  failure  mode  under  static  combined  uniaxial  compression 
and  shear  loading. 

7.  For  metal  panels  under  combined  compression  and  shear  loading 
stiffener  crippling  was  the  dominant  failure  mode.  Permanent 
deformation  of  the  skin  was  either  concurrent  or  precipitated  by 
stiffener  crippling, 

8.  Under  constant  amplitude  fatigue  loading  metal  panel  failure 
occurred  by  crack  initiation  in  the  skin  adjacent  to  a 
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stiffener.  The  static  strength  of  the  metal  panels  was 
unaffected  by  the  skin  crack  propagating  across  an  entire  bay. 
The  crack  initiation  life  was  approximately  60  percent  of  total 
number  of  cycles  sustained  by  the  panel  prior  to  skin  rupture  in 
a  single  bay.  The  nature  of  the  combined  loading,  i.e. 
compression  dominated  or  shear  dominated,  did  not  affect  the 
crack  propagation  pattern..  However,  the  initiation  sites  in  the 
two  cases  were  different.  Under  shear  dominated  loading  the 
cracks  initiated  at  stiffener  attachment  fastener  holes,  whereas 
under  compression  dominated  loading  the  cracks  initiated  in  the 
skin  at  the  edges  of  the  stiffener  flange  attached  to  the  skin. 

9.  Durability  c.  nsiderations  can  severely  limit  the  postbuckled 
operation  range  of  metal  panels.  In  the  panel  design  tested  the 
static  strength  range  was  250  percent  of  the  average  initial 
buckling  load.  However,  for  a  100,000  cycle  constant  amplitude 
fatigue  life,  the  panel  loads  would  have  to  be  restricted  to  125 
percent  of  the  initial  buckling  load. 

10.  Composite  panels  demonstrated  a  high  fatigue  threshold  relative 
to  the  initial  skin  buckling  loads.  Composite  panels  designed 
for  a  static  strength  equal  to  250  percent  of  the  initial  skin 
buckling  load  can  be  safely  operated  under  fatigue  loading  up  to 
200  percent  of  the  initial  buckling  load. 

11.  Composite  panels  tested  in  the  program  showed  a  greater 
sensitivity  to  shear  dominated  fatigue  loading  as  compared  with 
compression  dominated  fatigue  loading. 

12.  The  fatigue  failure  mode  in  composite  panels  was  separation 
between  the  cocured  stiffener  and  the  skin.  In  particular,  the 
region  at  Che  intersection  of  the  stiffener  and  the  ring  was 
vulnerable  to  the  failure  mode . 

13.  Repeated  buckling  had  no  influence  on  initial  skin  buckling 
loads  for  either  the  composite  or  the  metal  panels. 

14.  The  semi-empirical  design  methodology  was  used  to  develop  a 
design  procedure  for  composite  and  metal  panels  under  combined 
uniaxial  compression  and  shear  loading. 
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1.  A  single-mode  and  multi-mode  energy  method  based  postbuckling 
analysis  was  developed. 
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APPENDIX  A 
ANALYSIS  DETAILS 


This  Appendix  defines  the  individual  energy  integrals  used  in 
Equations  (54)  through  (57) ,  details  the  nonlinear  system  given  by  equation 
(58)  and  presents  closed  form  expressions  of  all  the  energy  integrals. 


The  displacement  functions  are  expressed  in  a  general  form  in 
Equation  (53) .  The  following  notations  are  used  for  the  derivatives  of  the 
displacement  functions: 
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The  individual  integrals  F,  G,  H  and  I  in  equation  (54)  are  defined 
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The  integrals  J  and  K  appeared  in  Equations  (55)  -  (57)  are  defined 
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The  system  of  nonlinear  algebraic  Equations  (58)  are  obtained  by 
setting  the  derivative  of  the  total  potential  energy  with  respect  to  each  of 
the  unknown  coefficient  to  zero.  Six  groups  of  equations  are  obtained  and 
they  are  given  below.  In  the  following  equations  the  total  load,  Pxx  and  PXy 
are  used  in  place  of  bNxx  and  bNXy. 
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r  .  *  00  ,  *  33  D*  34  43 

+  Gnm  ^  A22  G34runij  +  — ^  ^11  G34nmij  +  ~ (G34nmij  +  G34nmij ) 

2D*  35  53  >  *  44 

+  (G34nmij  +  G34nmij )  +  ~r“  ®22  G34nmij 

az  b-3 

2D*  45  54  4D*  55 

+  (G34nmij  +  G34nmij )  +  ~ G34nmij 


bw  *  011  101  110 

+  ~ a12  (^345runij  +  ^345mnij  +  ^345nmij ) 


aw  *  022  202  220 

+  r~  A22  (^345nmij  +  ^345nmij  +  ^345nmij ) 


3bw2  *  1111  w2  *  *  1122 

- 1  A11  I3455nraij  +  ~  <A12  +  2A66>  (I3455nmij 


1212  2112  22113  3aw2  2222 

+  2^3455nmij  +  2^-3455nmij  +  ^3455nmij)  +  g  A22  *3455nmij  | 

*  00  *  33  D*  34  34 

m~~~  a22  G44ijnm  +  — r  ^11  G44ijnm  +  — r~  (G44ijnni  +  G44nmij) 

D4  o-i  3.0 
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2D*  35  35  aD*  44  2D*  44 

+  _1£  (G44ijnm  +  G44nmij)  +  —2£  G44ijnm  + - ^  (G44ijnm 


45  4D*  55  bw  *  Oil  Oil 

+  G44nmij)  +  ~g~  G44ijnm  +  A12  ^445ijnm  +  ^445nraij 


110  aw  *  022  022  220 

+  ^445^,^)  +  — p  A22  (H445ijnm  +  ^445nmij  +  ^445ijnm) 


3bw2  *  1111  w2  *  *  1122 

+  - V  A11  I4455ijnm  +  “  <A12  +  2a66>  (^SSijnn 


1212  1212  2211  3 aw2  *  2222 

+  2l4455ijnm  +  2*4455nmij  +  ^4455ijnm)  +  °  A22  *4455ijnm 


+  ax 


11  aw  *  22 

G45ij  +  ~  a12  G45ij 


^  21  12 

+  bi  £  wo  a66  (G45ij  +  G45ij)  j 


r  .  *  111  ,  *  122 

"  *  Anm  Gpq  ~ ~  A11  ^134nmpqij  +  g  A12  ^134nmpqij 


*  212  221 
+  ~  a66  (Hl34nrapqij  +  ^134nmpqij) 


r  .  *  111  1  *  122 
'  Anm  Dpq  — —  A11  ^144nmijpq  +  £  A12  ^144nmijpq 


*  212  212 
+  —  a66  (^144nmijpq  +  ^144nmpqij) 


r  .  *  211  *  222 
-  Bjm,  Cpq  A12  H234runpqij  +  ^  A22  H234nmpqij 


*  112  121 

+  -  Agg  (H234nmpqii  +  ^234nmpqij) 

a 
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r  1  *  211  *  222 

®nm  Dpq  A^2  H244nmijpq  +  —  &22  ^244nmijpq 


,  *  112  112 
+  »  Ag6  H244nmijpq  +  H244nmpqij) 


r  h  *  011  110 

^nm  ^pq  A12  (2H334nnipqij  +  H334nmpq£j) 


*  022  220  3bw  *  1111 

+  rrr  A22  (2H334nmpq^j  +  H334rUBpqjj  )  +  — An  ^3345runpqij 


w  *  *  1122  1212 

+  (A12  +  2Agg)  (l3345nmpqij  +  2l3345nmpqij 


1221  2211  3 aw  *  2222 

+  2l3345nnipqj[j  +  l3345nmpqij)  +  - -  A22  l3345nmpqij 


r  *  011  110  110 

^nm  Dpq  A12  (H344nmijpq  +  H344nmijpq  +  H344Iunpqij  ) 

*  022  202  202 
+  - TjT  A22  (^344nmijpq  +  ^344mnijpq  +  ^344nmpqij) 


3bw  *  1111  v  *  *  1122 

— ~~  An  *3445mnijpq  +  ~r  (A12  +  ^Ags)  (l3445nmijpq 


1122  2112  1221  2211 
+  *3445runpqij  +  ^3445nmijpq  +  ^3445nmijpq  +  *3445nmijpq 


2211  3 aw  *  2222 

+  *3445nmpqij)  +  — A22  ^3445nmijpq 


r  .  *  Oil  Oil  *  022 

^run  ^pq  A12  (^444ijnmpq  +  2H444nmijpq)  +  A22  (H444ijnmpq 


022  3bw  *  1111 

+  2H444nmijDq)  +  - -S  An  l4445ijnmpq 

2aJ 
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w  *  *  1122  1122  1221 

+  2^  +  ^^4445ijtiinpq  +  *4445iimpqij  +  *4445ijranpq 


1221 


3aw 


*  2222 


+  2l4445mijpq)  +  A22  *4445ijnmpq  j 


&\  C 


rail 


*  11 


^  22 

An  ®34nmij  +  ~~  A12  ®34nraij  j 


[,  *  11  "  i.C  -1 

— -  An  ^44ijran  +  ~~  A12  ^44ijnm 


*  22 


r  *  21  21 
*>1  ^nm  Agg  (^34ijnm  +  ®34nmij) 


r  *  12  12 

M  Dnm  A66  (^44ijran  +  G44nmij) 


crim  c 


*  llll 


1122 


pq  crs  A11  ^ISAranpqrsij  +  7“  (Ai2  +  2A66)  (I3334ranpqrsij 


1.221  *  2222 

^  ^3334ranpqrsij ^  +  f~  A22  *3334nmpqrsij 

abJ  J 


*  1111 


1122 


••ran 


r  o-L.  JLJ-ii  ,  "  1  iii 

Cpq  Drs  An  l3344nmpqijrs  +  7^  (A12  +  2Agg)  (l3344nmpqijrs 


12i2  1212  2211 

+  2l3344ranpqijrs  +  2l3344tunpqrsij  +  ^3344ranpqijts) 


+  3a 
2a3 


*  2222 


A22  *3344ranpqijrs 


*  1111 


OL  XiXi. 

cran  ^pq  ®rs  ^  An  l3444mijpqrs 


*  *  1122  1122  2112 

+  — —  (A12  +  2Agg) (l3444ranijpqrS  +  2l3444nmpqijrs  +  2l3444nmijpqrs 
2ab 
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2112  -  *  2222 
+  i3444nnipqrsij )  +  22  *3444nmijpqrs  | 


r  h  *  1U1 

^nm  ^pq  ^rs  ~ J  A11  *4444ijnmpqrs 


+  -i- 

2ab 


*  *  1122  1122 
(A12  +  2Agg) (l4444ijnmpqrs  +  *4444nmpqijrs) 


+ 

2b3 


*  2222 

22  ^-44441 jnmpqrs 


r  ,  *  00  ,  *  33  D*  34  43 

*  W°  [  ^2  A22  G45i:J  +  3  Dl1  G45i3  +_ab  (G45iJ  +  G45iJ) 

2D*  35  53  *  44  2D*  45  54 

+  (G45ij  +  G45ij)  +  -S-  D22  G45ij  +  — (G45ij  +  G45ij) 

bJ  b^ 


+  it  %6  G««  ] 

2  r  .  *  011  110 

w°  I  iSn  Al2  (H455iJ  +  2H455ij) 


3  r  *  1111  *  *  1122  2112 

'  W°  L  I3  Al1  145551:5  +  2lb  <Al2  +  2Afi6)  (l4555ij  +  I^555ij) 


*  2222 

+  A22  l4555ij 

ab3  J 


_5H_  =  0 

da.\ 


r  *  1  1 

Anm  b  Ml  Mnm  +  As^s  Mnm 


it  2 

+  ®mn  [  a  M2  ^2nm 
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+  Gnm 

r  ,  *  0  bw  *  11  aw  *  22  n 

““  Al2  F3nn)  +  _o_  All  G35nm  +  Al2  G35nm  1 

+  bnm 

r  K  *  0  bw  *  11  *  aw  *  22  n 

~~  a12  F4nm  +  A11  G45nm  +  A12  G45nm 

* 

+  ai  (ab  An  +  a  AsEs) 

"  Fxx  s 

1  *  Grmi  cpq 

r  .  *  n  *  22 

~2^  ^33nmpq  +  ^  A12  G33nrapq 

'  Gnm  bpq 

•k  ^  22 

-fe-  An  G34nmpq  +  A12  G34nmpq  j 

*  11  #  22 
Dnm  Dpq  f  ~~  A11  G44nmpq  +  ~~  A12  G44runpq 


abw  *  0  bw^  *  11  aw^  *  22 

-  a12  f5  -  -jf  A11  g55  -  a12  g55 


-22-  _  o 

9b! 


*  2 

Anm  (a  a66  Flnm) 


*  1 

+  Enm  (b  a66  F2nm> 


r  *  21  12 

+  Gnm  wo  a66  (G35nm  +  G35nnr 


r  *  12  21 

+  bnm  wo  a66  (G45run  +  G45nnr 
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*  r  12  12  21 

a  ?xy  "  ^66  ^nm  ^pq  ^33nmpq  +  ^nm  ^pq  ^34nmpq  +  ^34nmpq) 


12  12 
+  ^nm  Dpq  G44nmpq  +  w0  G55 
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The  individual  energy  integrals  defined  in  Equations  (A- 2)  through 
(A- 8)  are  evaluated  in  closed  form  for  the  particular  displacement  functions 
given  by  Equations  (59).  These  integrals  are  expressed  in  terms  of  simple 
integrals  of  the  sine  and  cosine  functions  and  the  combinations  of  those 
functions.  The  assumed  displacement  functions  are  given  in  equations  (59)  in 
Section  3  and  they  are  rewritten  below. 

fl(n,m;x,y)  -  X2  (l,n)  cos  mwy  +  X^  (l,n)  sin  mwy 

f2(n,m;x,y)  -  [cos  nwx  -  X3(l,n)]  cos  nmy  +  [sin  nwx  -  X2(n,l)]  sin  mwy 

f3(n,m;x,y)  -  sin  nwx  sin  mwy 
f4(n,m,x,y)  -  Xi(l,n)  Y2(l,m)  -  X2(l,n)  Yi(l,m) 
f5(x,y)  -  sin  wx  sin  wy 
where 

Xi(n,m)  -  sin  nwx  sin  mwx,  Y^(n,m)  -  sin  nwy  sin  mwy 

X2(ntm)  -  sin  nwx  cos  mwx,  Y2(n,m)  -  sin  nwy  cos  mwy 

X3(n,m)  -  cos  nwx  cos  mwx,  Y3(n,m)  -  cos  wry  cos  mwy 

The  derivatives  of  the  displacement  functions  are  given  by 

fl,l  (n,m;x,y)  -  w[Zi(n)cos  mwy  +  Z2(n)  sin  nwry] 

f^,2  (n,m;x,y)  -  mw[Xi(l,n)cos  mwy  -  X2(l,n)  sin  rowy] 

f2,i  (n,m;x,y)  -  w{[n  X2(n,l)  +  X2(l,n)  -  n  sin  nwxj  cos  miry 

+  [n  cos  nwx  -  n  X3(l,n)  +  Xi(l,n)]  sin  mwy} 


A 


4 


f2,2  (n,m;x,y)  =  mw{[sin  nwx  -  X2(n,l)]  cos  mwy 

-  [cos  nwx  -  X3(l,n)]  sin  mwy)  v 

k 
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f 3 > l  (n,m*x,y)  -  nir  cos  nirx  sin  miry 

f3,2  (n,m;x,y)  -  mir  sin  nirx  cos  miry 

f3,3  (n,m;x,y)  -  -n2ir2  sin  nirx  sin  miry 

^3.4  (n,m;x,y)  -  -m2ir2  sin  nirx  sin  miry 

£3,5  (n,m;x,y)  -  nmir2  cos  nirx  cos  miry 

f4-l  (n, m;x,y)  -  ir[Z2(n)Y2(l,m)  -  Zi(n)  Yi(l,m)] 

f4-2  (n,m;x,y)  -  ir[Xi(l,n)Zi(m)  -  X2(l,n)  Z2(m)] 

f4,3  (n,m;x,y)  -  n2{ [2n  X3<l,n) -  (1+n2)  Xi(l,n)]  Y2  (l,m) 

+  [2n  X2(n,  1)  +  (1+n2)  X2(l,n)J  Yxd.m)} 

f4>4  (n,m;x,y)  -  -  ir2{Xi(l,n)  [ (l+m2)Y2(l,m)  +  2m  Y2(m,l)] 

+  X2(l,n)[2m  Y3(l,m)  -  (1+m2)  Yi  (l,m)]) 

f4,5  (n,m;x,y)  -  ir2[Z2(n)Z1(m)  -  Z*(n)  Z2?m)] 
f5>l  (x,y)  -  7r  cos  7rx  sin  iry 

£5, 2  (x,y)  -  ff  sin  irx  cos  iry 

f5 ,3  (x,y)  -  -it2  sin  irx  sin  iry 
£5-4  (x,y)  -  -7r2  sin  irx  sin  iry 
f5'5  (x,y)  "  ir2  cos  irx  cos  iry 
where , 

ZX  (n)  -  X3  (l,n)  -  n  Xj.  (l,n)  ,  \  (n)  -  Y3  (l,n)  -  n  Yx  (l.n) 

Z2  (n)  -  X2  (n, 1)  +  n  X2  (l,n)  ,  Z2  (n)  -  Y2  (n.l)  +  n  Y2  (l.n) 

The  integrals  are  expressed  in  terms  of  a  series  of  elementary  integrals  and 
combination  of  these  elementary  integrals  defined  below. 

1  0  for  n  *  0 

I(n)  •*=  cos  nirxdx  -  • 

,0  l  1  for  n  -  0 
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J 


Io(n,m)  - 


Jl(n,m)  - 


l2(n,m)  - 


Jo(n,m,p) 


^(n.rn.p) 


JaCn.m.p) 


l3(n,m,p) 


(n) 


1  0  for  even  n 


sin  njrxdx  -  •< 


-2—  for  odd  n 
nn 


0 

,1 

0 

,1 


sin  njrx  sin  nwrxdx  -  -i-  [I(n-m)-I(n+m)  ] 

2 


sin  rorx  cos  mjrxdx  -  — [J(n-m) -J(n+m)  ] 


cos  njrx  cos  mjrxdx  -  — ±—  [I(n-m)+I(n+m) ] 

2 


sin  njrx  sin  nwrx  sin  pjrxdx 


-i-  [ J (n-m+p)+J (n+m-p) -J (n-m-p) -J (n+m+p) ] 

4 


sin  njrx  sin  mjrx  cos  pjrxdx 


-i-  [I(n-m+p) -I(n+m-p)+I (n-m-p) -I(n+m+p) ] 
4 


sin  njrx  cos  mjrx  cos  pjrxdx 


-i~  [  J (n-m-!-p)+J (n+m-p)+J  (n-m-p)+J  (n+m+p)  ] 
4 


cos  njrx  cos  mjrx  cos  p7rxdx 


-i-  [I(n-m+p)+I(n+m-p)+I(n-m-p)+I(n+m+p) ] 
4 
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Io(n,ni,p,q)  - 


Jl(n,m,p,q)  - 


l2(n,m,p,q)  - 


^(n.m.p.q)  - 


l4(n,m,p,q)  - 


Jo(n,m,p,q,r) 


fl 

sin  nrrx  sin  mjrx  sin  pwx  sin  qjrxdx 

.  0 

-  [Ii(n,m,p-q)  -  Ii(n,m,p+q)] 

fl 

sin  n?rx  sin  nwrx  sin  p?rx  cos  qwxdx 

.  0 

—  [Jo(n.m,p-q)  +  Jo(n,m,p+q)] 

fl 

sin  njrx  sin  mjrx  cos  prrx  cos  ^/rxdx 

.  0 

-i-  [Ii(n,m,p-q)  +  Ii(n,m,p+q)] 

2 

fl 

sin  nirx  cos  mwx  cos  pjrx  cos  qjrxdx 

.  0 

~  [J2(n,m,p-q)  +  J2(n,m,p+q)] 

fl 

cos  n7rx  cos  nwrx  cos  pjrx  cos  qwxdx 

.  0 

[l3(n,m,p-q)  +  l3(n,m,p+q) ] 

1 

-  sin  wrx  sin  nmx  sin  p7rx  sin  qwx  sin  rjrxdx 
.  0 

- -  [Jl(n,m,p,q-r)  -  Ji(n,ni,p ,q+r)  ] 
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Il(n,m,p 


J2(n,m,p 


l3(n,ra,p 


J4(n,m,p 


l5(n,m,p 


Io(n,m,p 


-q,r) 


sin  nrrx  sin  nwrx  sin  prrx  sin  qrrx  cos  r^xdx 


y-  [Io(n,m,p,q-r)  +  Io(n,ra,p,q+r) ] 


.q,r) 


sin  n?rx  sin  nwrx  sin  p7rx  cos  qrcx  cos  rrrxdx 

0 

-i-  [Ji(n,m,p,q-r)  +  Jx(n,m,p,q+r) ] 

2 


>q.r) 


,q,r) 


sin  nxx  sin  nwrx  cos  prrx  cos  qffx  cos  rwxdx 


-i-  [l2(n,m,p,q-r)  +  l2(n,m,p,q+r) ] 


sin  nrrx  cos  nwrx  cos  prrx  cos  qrrx  cos  rjrxdx 


~~  [J3(n,m,p)q-r)  +  J3(n,m,p,q+r) ] 


,q,r) 


cos  nwx  cos  nwrx  cos  p7rx  cos  qwx  cos  rrrxdx 


0 


- i-  [l4(n,m,p,q-r)  +  l4(n,m,p,q+r) ] 


q,r,s) 


fl 

sin  n?rx  sin  nwrx  sin  prrx  sin  qrrx  sin  rrrx  sin  srrxdx 

.  0 

-L.  [I1(n,m,p,q,r-s)  -  Ix(n,m,p,q,r+s) ] 
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JlCn.m.p.q.r.s)  - 


sin  n?rx  sin  nwrx  sin  pxx  sin  qjrx  sin  rjrx  cos  sjrxdx 


-  ~  (Jo(n,m,p,q,r-s)  +  Jo(n,ra,p,q,r+s) ] 


l2Cn>m,P1q»r.s) 


sin  njrx  sin  nwrx  sin  pjrx  sin  qjrx  cos  rjrx  cos  sjrxdx 


- [Ii(n,m,p,qlr-s)  +  IiCn.m.p.q.r+s) ] 


J3(n,m,p,q,r,s)  - 


sin  nrrx  sin  nwrx  sin  pjrx  cos  qjrx  cos  r?rx  cos  sjrxdx 

0 

[J2(ti,ra,p,q,r-s)  +  J2(n,m,p,q,r+s) ] 


,1 


I^Cn.m.p.q.r.s)  - 


sin  njrx  sin  nwrx  cos  pjrx  cos  qjrx  cos  rjrx  cos  sjrxdx 

0 

y-  [^(n.m.p.q.r-s)  +  I3(n,m,p,q,r+s) ] 


JsCn.m.p.q.r.s) 


sin  n?rx  cos  nurx  cos  p7rx  cos  q7rx  cos  rjrx  cos  sjrxdx 


l6(n,m,p,q,r,s) 


- [J4(n,m,p,q,r-s)  +  ^(n.m.p.q.r+s) ] 


cos  njrx  cos  nwrx  cos  pjrx  cos  qjrx  cos  rjrx  cos  sjrxdx 


—  [l5(n,m,p,q,r-s)  +  l5(n,m,p,q,r+s) ] 
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In  addition  to  the  elementary  integrals  defined  above,  the  following 
combinations  of  elementary  integrals  are  also  defined. 

f1 

Bi(i,j,k,0  -  cos  iwx  cos  jwx  Z2(k)  Z2(Odx 

Jo 

-  I4(k,t,i,j,l,l)  +  kl4(l,t,i,j,k,l) 

+tl4(l,k,i, j ,4,1)  +  ktl4(l,l,i,j,k,t) 

fl 

B2(i,j,k,t)  -  cos  ijrx  cos  jwx  Zi(k)  Z2(Odx 

Jo 

-  J5(t,i,j,l,l)  -  kJ3(l,k, t , i.j ,1) 

+iJ5(l,i,j ,k,i,l)  -  kiJ3(l,l,k,i,j , t) 

fl 

B3(i,j,k,t)  -  I  cos  ijrx  cos  jirx  Z^(k)  Zi(i)dx 

Jo 

-  I6(i,j,k,t,l,l>  -  kl4(l,k,i,j,i,l) 

-tl4(l,*,i,j,k,l)  +  ktl2(l,l»k, t ,i, j) 

fl 

B4(i,j,k,t)  -  sin  iirx  sin  jirx  Z^(k)  Zi(t)dx 

Jo 

-  I4(i. j ,k, t ,1,1)  -  kl2(l,i, j ,k, t ,1) 

-tl2(l,i,j,t,k,l)  +  ktIo(l,l,i,j,k,t) 

f1 

B5(i,j,k,t)  -  I  sin  iwx  sin  j?rx  Z^(k)  Z2(Odx 

Jo 

=  J3(i,j,t,k,l,l)  -  kJid.i.j.k.t.l) 

+tJ3<l,i,j ,k, l  ,1)  -  ktJi(l,l,i,j,k, t) 
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B6(i,j,k,0 


Ci(i,j,k,0 


C2(i,j ,k,i) 


C3(i,j.k,0 


C4(i.j.k,0 


Did.j.k,*) 


D2(i,j ,k,t) 


-  sin  iirx  sin  jirx  Z2(k)  Z2(t)dx 
,  0 

-  I2(i. j .k, t ,1,1)  +  kl2(l,i,j ,4,k,l) 
+4l2(l,i,j,k,t,l)  +  ktl2(l,l,i, j ,k, t) 

f1 

-  I  sin  iirx  cos  jirx  Xi(l,k)  Z2(t)dx 

Jo 

-  I2(l»i.k, t , j ,1)  +  il2(l,l.i,kfjt*) 

fl 

-  sin  iirx  cos  jirx  Xi(l,k)  Z]_(i)dx 
.  0 

-  J3(l,ifk,jIt,,l)  -  tJid.l.i.k.t.j) 

r1 

-  I  sin  iirx  cos  jirx  X2(l,k)  Z2<t)dx 

Jo 

-  Jsd.i.t.j.k.l)  +  tJ3(l,i,k,j,t,l) 

f1 

-  j  sin  iirx  cos  jirx  X2(l,k)  Z^(t)dx 

Jo 

-  I4<l,i,j,k,t,l)  -  tl2(l,l,i,t,j ,k) 

fl 

-  sin  iirx  X2(l,j)  X2(l,k)  X2(l,t)dx 
.  0 

-  _1_  [l3(l,i»j .k, O  -  I4(l, i, j ,k, 4,2)] 

rl 

-  sin  iirx  X^(l,j)  X2(l,k)  X2(l,t)dx 

Jo 

-  -j-  [J2(l.i.j.k,0  -  J3(l.i.  j  .k,  t  ,2)  ] 
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D3(i,j,k,0  -  I  sin  i*x  Xjd.j)  X1(l.k)  X2(l,t)dx 

Jo 

-_1_  [Iid.i.j.k.i)  -  I2(l»i»j  »k,  t  ,2)  ] 

D4d>j.k,t)  -  sin  ifrx  X^d.j)  X^d.k)  Xi(l,t)dx 

Jo 

*~  ^  -  [J()d >  ^ >  j  >k,  t )  -  J]_(l,  ii  j  fk,  t  ,2)  ] 

r1 

Ei(i,j,k,t)  -  j  cos  iirx  Z2<j)  Z2(k)  Z2(t)dx 

Jo 

-  -i-  {J2(j,k,t,i,l)  +  J3<j,k>t,i,l,2)  +  j  J3(2 ,k, t , i , j , 1) 

+  k  J3(2,j,t,i»k,l)  +  :  J3(2,j ,k,i,t,l)  +  jk  J3d,2,t ,i, j ,k) 
+  jt  J3(l,2,k,i»j  » O  +  kt  J3(l,2,j ,i.k,t) 

+  jkt  {J4(l)i,j,k,0  -  J5d,i,j.k,t,2)J} 

E2<i.j.k,t)  -  I  cos  iirx  Zi(j)  Z2(k)  Z2(t)dx 

Jo 

-  d3(k,t,i,j,l)  +  I4(k, t ,i,j ,1»2)  -  j  I2(2,j,k,t,i,l) 

+  k  I4(2, 4 ,i, j ,k,l)  +  «•  I4(2,k,i,j,*,l)  -  jk  I2(l,2,j , t ,i>k) 

-  jt  I2d,2.j,k,i.O  +  kt  I4(l,2,i»  j  »k,  t) 

*  jkt  [ I3 (1 » j .i.k.t)  -  I4(l, j ,i»k, t ,2) ] } 

E3(i,j,k,t)  -  I  cos  i7rx  Zx(j)  Zi(k)  Z2(t)dx 

Jo 

-  -i-  {J4(t,i.j  «k,l)  +  J5(t,i,j,k,l,2)  -  j  J3(2,j , t ,i,k,l) 

-  k  J3(2 ,k, t , i» j  » 1)  +  <■  J5(2,i,j,k,t,l)  +  jk  Ji(l, 2 , j ,k, t , i, ) 

-  jt  J3(l,2,j ,i,k,  t)  -  kt  J3(l,2,k,i,j . O 
+  jkt  [J2(l,j  ,k,i.  i)  -  J3d»j»k,i,t,2)]} 
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E4(i,j.k.O 


Fi(i,j,k,t) 


F2(i,j,k,0 


F3(i,j ,k, i) 


F4(i,j ,k, i) 


F5(i,j ,k,t) 


I  cos  i irx  Zx(j)  Zj.(k)  Z^ ( 4 ) dx 
0 

-i-  {I5(i, j ,k, 4,1)  +  I$(i, j ,k, 4 ,1, 2)  -  j  I4(2, j , i.k, t ,1) 
k  I4(2,k, i, j ,4.1)  -  i  I4(2, t ,i, j ,k,l)  +  jk  I2(l,2, j ,k,i, O 

jt  I2(1|2,  j  ,  i  +  k*  I2(l|2  ,k,  4  ,  i.  j ) 

jkt  [Ii(l,j,k,4,i)  *  I2(l. j  »k, * , i,2) ] } 


(cos  i  irx  Z2(j)  X3(l,k)  X^(l,4)dx 
0 

-L-  {J!(l,2, j ,k, t , i)  +  j[J2(l,k,4,i,j) 

fl 

cos  iirx  Z2(j)  Xi(l,k)  X2(l,4)dx 

.  0 

-i-  {I2(l,2,j  ,k,i,  4)  +  j[I3(l,k,i,j,0 

fl 

cos  iirx  Z2(j)  X2(l,k)  X2(l,4)dx 

.  0 

-L.  { J3  (1 , 2 ,  j  ,i,k,  4)  +  j[J4(l,i,j,k,4) 

r1 

cos  iirx  Z^(j)  X3(l,k)  X^(l,4)dx 

.  0 

-i-  {I2(l,2,k, 4 ,i. j)  -  j[Ii(l,j,k,4,i) 


J3(l,k,4,i,j,2)]} 


14(1. k,  i,j  .  *•  .2)  ] } 


J5(l,i,j,k,4,2)]} 


I2(l,j,k,4,i,2)]} 


.1 

cos  iirx  Z^(j)  X^Cl.k)  X2(l,4)dx 

.  0 

-i_  {J3(l,2,k,i.j . O  -  j [J2(l,j ,k,i, 4)  -  J 3 ( 1 > j ,  k ,  i ,  4 , 2 ) ] } 
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F6(i,j,k, 


Gid.j.k, 


G2(i.j ,k, 


G3(i, j ,k, 


G4 ( i , j ,k,t 


-  cos  iwx  Zx(j)  X2(l,k)  X2(l,t)dx 

JO 

-  -i-  {14(1,2,1, j  >k, t)  -  j(I3(l,j,i,k,0  -  I4(l,j,i, 

r1 

-  sin  Ijtx  X2<l,j)  Z^(k)  Z^(t)dx 

Jo 

-  -i-  {14(2,1, j,k,t,l)  -  kl2(l,2,i,k,j ,t)  -  tl2(l,2, 

+  ktdid.i.k.t.j)  -  I2(l,i,k,  t , j  ,2) ] } 

-  I  sin  iirx  X2(l,j)  Z^(k)  Z2d)dx 

Jo 

-  -L-  {J3(2 , 1 , t , j ,k, 1)  -  kJ1(l,2,i,k,t,j)  +  tJ3(l,2, 

-  kt[J2(l,i,k,j ,t)  -  J3(l,i,k, j ,4,2)]} 

fl 

sin  ijrx  X2(l,j)  Z2(k)  Z2d)dx 

Jo 

-  _i_  {I2(2,i,k, t , j ,1)  +  kl2(l,2,i, 1 ,j ,k)  +  4l2(l,2,J 

+  ki[I3(l,iJ.k.O  -  I4(l,i,  j  ,k,  t  ,2)  ] } 

rl 

-  sin  i?rx  Xi(l,j)  Z^(k)  Z^(i)dx 

Jo 

-  -i-  {J3(2, i, j ,k, 4,1)  -  kJi(l,2,i,j,k,i)  -  tJi(l,2,i 

+  kt[Jo(l,i,j,k,t)  -  Ji(l,i,j,k,t,2)]} 


,4,2)]} 


,4,j,k) 


.j.k.O 


,k, j , t) 


,j  ,4, k) 
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Gsd.j.k.O  -  sin  i*rx  Xi(l,j)  Zi(k)  Z2d)dx 
JO 

- -i-  {I2(2,i,j,t,k)l)  -  kI0(l,2,i,j,k,O  +  tl2(l,2,i,j,k,0 

-  kttlxd.i.j.k.O  -  I2(l,i,j  »k,t,2)]} 

G6(i,j,k,t)  -  sin  ixxXid.j)  Z2(k)  Z2(t)dx 

Jo 

-  -i-  lJi(2,i,j,k,i,l)  +  kJxd^.i.j.t.k)  +  *Jid,2,i,J.k.O 

+  kt [J2(l,i. j  ik, t)  -  J3d,i,j,k>t,2)]} 

Ko(i,j,k,0  -  Z2(i)  Z2(j)  Z2(k)  Z2(Odx 

Jo 

-  -I_  {3I0(i,  j  ,k,  4  , )  +  4l1(i,j,k,l,2)  +  ^(i.j.k.t.A) 

8 

+  i[2Ii(2,j,k,t,i)  +  Ii(4J,k,t,i)]  +  jddd.i.k.i.j) 

+  Ii(4,i,k,t,j)]  +  kdixdd.j.t.k)  +  I1(4,i,j,t,k)] 

+  t  [  2Ii  (2 ,  i ,  j  ,k,  i)  +  Ii(4d,j,k,t)]  +  ij[I2(ktl[i,j) 

-  I3(k,t>i,j,4)]  +  ik[I2(j , t ,i,k)  -  I3(j ,  t , i»k,4) ] 

+  i*[I2(j,k,i,0  -  I3(j  »k,i»  i  >4)  ]  +  jk[I2(i,t>j,k) 

-  I3(i»*»j.k,4)]  +  jt[I2(i,k,j,0  -  I3(i,k,j , t ,4)] 

+  kt[I2(isj ,k,t)  -  I3(i, j ,k, t ,4) ]  +  ijk[2I3(2, i ,i. j ,k) 

-  I3(4,  i ,  ii  j  .k)  ]  +  ijt[2I3(2,k,i,j,0  -  I3(4,k,iJ,Ol 

+  iki [2I3(2,j ,i>k, t)  -  I3(4, j ,i,k, O ]  +  jkt [2I3(2, i, j ,k, i) 

-  I3(4,ii  j  »k,  0  3  +  ijkt[3I4(i,j,k,0  -  4I5(i,  j  ,k,  i  ,2) 

+  l5(isj,k,*,4)]} 
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Ki<i,j,k,0  -  Z!<i)  Z2(j)  Z2(k)  Z2 ( O dx 

Jo 

-  -I-  OJiCj.k.t.i)  +  4J2(j ,k, 4,1,2)  +  J2(j,k,4,i,4) 

8 

-  i  { 2  J  o  ( 2 ,  i ,  j  ,k,  i)  +  J0(4,i,j,k,O)  +  j  [2J2(2,k,  4  ,i,  j) 

+  J2(4,k, t ,i,j)]  +  k[2J2(2, j , 4 ,i,k)  +  J2(4, j , 4 , l,k) ] 

+  t [2J2(2, j ,k, i, t)  +  J2(4,jfk,i,-)]  -  ij[Jl(i,k,t,j) 

-  J2(i,k,  t ,  j  ,4)  ]  -  iktJxU.j.i.k)  -  J2(i,  j  .  t  ,k,4)  ] 

-  i4[Ji.(i,j,k,t)  -  J2(i,j,k,4,4)]  +  jk[J3(4,l,j,k) 

-  J4(‘ »ii j  >k,4) ]  +  ji[J3(kfl,j,4)  -  J4(k,i,j,t,4)] 

+  k4[J3(j,i,k,4)  -  J4(j  ,i|k,  t  ,4) )  -  ijk{2J2(2,i,4,j  ,k) 

-  J2(4, 1, t , j ,k) ]  -  ijt[2J2(2,i,k,j ,t)  -  J2(4, l,k, j , t) ] 

-  ikt [2J2(2,1, j ,k, t)  -  J2(4, i , j ,k, i) ]  +  jkt [2J4(2 , i , j ,k, t ) 

-  J4(4,i.j,k,t))  -  ijkt [3J3(i, j ,k, t)  -  4J4(i,J ,k, t ,2) 

+  J4<i>j.k,i,4)]} 

K2(i,j,k,t)  -  Z!(i)  ZX(j )  Z2(k)  Z2( 4 )dx 

Jo 

-  _1_  {3l2(k,4,i,j.)  +  4I3(k, 4 ,i, j ,2)  +  I3(k, 4 ,i,j ,4) 

8 

-  i[2I1(2,i,kt  4 ,J)  +  I1(4,i,k,4,j)]  -  j [2Ii(2, j ,k,  4 ,i) 

+  I1(4,j,k,4,i)]  +  k[2I3(2,4,i,j ,k)  +  I3(4, 4 ,i, j ,k) ] 

+  4 [2I3(2,k, i, j , 4)  +  I3(4,k,i,j,4)]+  ij[I0(i,j,k,4) 

-  Ii(i, j ,k, 4,4)]  -  ik[I2(i,4,j,k)  -  I3(i»  4 , j ,k,4) ] 

-  i4[I2(i,k,j,4)  -  I3(i,k,j,4,4)]  -  jk[I2(J,4.i,k) 

-  I3(j , 4 ,i .k,4) ]  -  j4[I2(j,k,i,4)  -  I3(j ,k, i, 4 ,4) ] 

+  k4[I4(i,j,k,4)  -  I5(i,j,k,4,4)]  +  ijk[2I1(2,i,j,4,k) 

-  Ii(4,i,j,4,k)]  +  ij4[2I1(2,i,j,k,4)  -  I1(4,i,j,k,4)l 

-  ik4[2l3(2,i,k,j,4)  -  I3(4,i,k,j,4)}  -  jk4[2I3(2,j,i,k,4) 
.  I3(4,j ,i,k, 4)]  +  ijk4[3I2(i,j,k,4)  -  4I3(i, j ,k, 4 ,2) 

+  I3(i* j ,k, 4 ,4)] } 
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K3(i,j,k,t)  -  Z!(i)  Zi(j)  ZX(k)  Z2<Odx 

Jo 

-  JL  {3J3(t ,  i,  j  ,k)  +  4J4(t,i,j,k,2)  +  J4(t  ,1 ,  j  .k,4) 

8 

-  i{2J2(2,i, i , j ,k)  +  J2(4,i,t,j,k)]  -  j [2J2(2, j , t ,i,k) 

+  J2(4,j,t,i,k)]  ‘  k{2J2(2,k,  t  ,i,  j)  +  J2(4,k,t,l,j)] 

+  t.  [ 2J4<2 , i ,  j  ,k,  O  +  J4(4, i,  j  ,k, t) ]  +  ij  Pld,  j  ,  ‘  .k) 

-  J2<iJ.‘»k,4)]  +  ik[Ji(i,k,t,j)  -  J2(i,k,  t  ,j  ,4)  ] 

-  it [J3(i,j ,k, t)  -  J4(ifj.k,t,4)]  +  jk[Ji(j ,k, t , i) 

-  J2(j (k, t ,1,4) ]  -  jt[J3(j,irk,0  -  J4(j,i,k,t,4)] 

-  kt [J3(k,i, j , t)  -  J4(k,i,j,*,4)]  -  ijk[2J0(2,i,j ,k, i) 

-  Jo<4,i,j,k,0]  +  ijt[2J2(2,i,j,k,t)  -  J2(4,  i,  j  ,k,  t )  ] 

+  ikt [2J2(2,i,k, j , O  *  J2(4,i,k, j , t) )  +  jkt [2J2(2, j ,k, i, O 
+  J2<^.3.k,i,01-  ijktiaJid.j.k,*)  •  4J2(i,  j  ,k,  i  ,2) 

+  J2(i, j ik, t ,4) ] } 

K4(i, j ,k, O  -  ZL(i)  Zi(j)  Zx(k)  Zi(t)dx 

Jo 

-  -L.  {3I4(i, j ,k, t)  +  4I5(i,j ,k,t ,2)  +  I5(i, j ,k, i ,4) 

3 

-  i[2I3(2, i, j ,k, t)  +  I3(4,i, j ,k, i)]  -  j [2I3(2, j ,i,k, t) 

+  13(4, j , i,k, t) ]  -  k[2I3(2,k,i, j , t)  +  I3(4,k, i, j , O ] 

-  t[2I3(2,t,i,j ,k)  +  I3(4, l  ,i, j ,k)]  +  ij [l2(i. j .k, t) 

-  I3(i, j ,k, t ,4) ]  +  ik[I2(i,k,j,*>  -  I3(i,k, j , t ,4) ] 

+  it[I2(i,t,j,k)  -  I3(ii * . j  »k,4) ]  +  jk[I2(j,k,i,t) 

-  I3(j ,k, i, t ,4) ]  +  jt[I2(j,t,i,k)  -  I3(j,t,i,k,4)] 

+  kt[I2(k,i,i,j)  -  l3(k,i.i,j,4)]  -  ijk[2Ix(2 , i, j ,k, t) 

-  I1(4,i,j,k,t)]  -  ijt[2I1(2,i1j>t,k)  -  Ii(4,i.j,‘.k)] 

-  ikt[2I1(2,i,k,t,j)  -  Ii(4,i,k,t,j)]  -  jkt [2Ix(2, j ,k, 1 , i) 

-  Ii(4, j ,k, t , i) ]  +  ijkt [3I0(i. j ,k, O  -  4I1(i, j ,k, t ,2) 

+  Il(i,j ,k,t,4)J } 
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LoU.j.k.O  - 


Li<i,j,k.O  - 


L2d.j  ,k.O 


L3(i,j,k,0  * 


L4(i,j,k,t)  - 


MooU.j.k.O  “ 


X2(l,i)  X2d,j)  X2(l,k)  X2(l, t)dx 

Jo 

4-  [3I4(i,j,k,t)  -  4I5(i,  j  fk,  t  ,2)  +  I5(i,j,k, 4.4)] 

8 

Xi(l.i)  x2(l,j)  X2(X,k)  X2d,*)dx 

Jo 

4“  [3J3(iJ.k.4>  -  4J4(i.j,k,t,2)  +  J4(i,j,k,t,4)] 
o 

XXd,i>  XXd.j)  X2(l,k)  X2(l, t)dx 

Jo 

4-  [3l2d,j,k,t)  -  4I3(i,j  ,k,  t  ,2)  +  I3d,  j  ,k,  t  ,4)  ] 

8 

xxd,i)  Xid.j)  xxd,k)  X2(l, t)dx 

Jo 

4-  [3JXd,jj  »k,  <.)  -  4J2(i,  j  ,k,  t  ,2)  +  J2(i,4,k,t,4)] 

8 

.i 

XXd, i)  xx(l,j)  Xx(l,k)  XX(l,Odx 

0 

4-  [3l0(i,j,k,O  -  4Ix(i,  j  ,k,  1 ,2)  +  Ix(i ,  j  ,k,  t  .4)  ] 

O 

Z2(i)  Z2( j )  X2(l,k)  X2(l,Odx 

Jo 

4-  {l2(i,j,k,0  -  IX(i,  j  ,k,  1 ,4)  +  i[2I3(2,  j  ,i,k,  t) 

-  I3(4,  j  ,  i.k,  t)  ]  +  j[2I3(2,itj,k,0  -  I3(4,i,j  ,k,  O] 
+  ij[3I4(i,j,k,t)  -  4I5(i, j ,k, 1 ,2)  +  Is(i, j ,k, t ,4) ] } 
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,k, O 


Z2(i)  Z2(j)  XxU.k)  X2(l,t)dx 


MioU.j  .k*0 


Mn(i,j,k,0 


Mi2(i,j,k,0 


-  4"  {Jl(i.j.k,0  -  J2(i,j,k,t>4)  +  i[2J2(2, j  ,k,  i,  t) 

o 

-  J2(4,j,k,i.O]  +  j[2J2(2,i,k,j,t)  -  J2(4,i.k,j  ,  t)  ] 
+  ij [3J3(k,i, j ,4)  -  4J4(k,i,j,*,2)  +  J4(k,i,j,t,4))} 

-  Z2(i)  Z2(j)  Xxd.k)  Xi(l,Odx 
JO 

--I"  do(i.j.k,t)  -  Ii(i,j,k,t,4)  +  il2Ii(2,J,k,i,i) 

O 

-  Ii(4, j ,k, t ,i) ]  +  j[2I1(2,i,k,t,j)  -  I1(4,i»k,t,j)] 
+  ij  [3I2(k,  t , i, j  )  -  4I3(k,t,i,j,2)  +  I3(k,i,iJ,4)]} 

-  Zxd)  Z2(j)  X2(l,k)  X2(l, t)dx 
.  0 

“4"  <J3(j.i>k,0  -  J4(j  ,  i»k,  t  ,4)  *  i[2J2(2 ,  i,  j  ,k.  1) 

O 

*  J2<4,i,j,k,0]  +  j[241(2,i,j,k,t)  -  J4(4,i,j,k,0] 

-  ij[3J3d,j,k,0  -  4J4(i,  j  ,k,  t  ,2)  +  J4(i,J  ,k,*  ,4)] } 
fl 

Zxd)  Z2 ( j )  Xx(l,k)  X2(l,  t)dx 

.  0 

“4"  {l2(j.k,i,0  -  I3(j  »k,i,  t  ,4)  -  i^Ixd.i.j.k.O 

O 

-  Ii(4,i,j,k,0]  +  j[2I3(2,k,i,j,0  -  I3(4,k,i.J,0] 

-  ij [3I2(i,k, j , 1)  -  4I3(i,k, j , t ,2)  +  I3(i,k, j , * ,4) ] } 

Z]d)  Z2(j )  Xxd.k)  Xxd.Odx 

Jo 

4-  {Jxd.k.td)  -  J2(j,k,t,i,4)  -  i[2J0(2,i,j,k,O 

O 

-  J0(4, i. j ,k, t) ]  +  j[2J2(2,k,t,i,j)  -  J2(4,k,i,i,j)] 

-  ij [3J1(i,k, t , j )  -  4J2(itk,*,j,2)  +  J2(i,k,t,j,4)J) 
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M20(i,j,k,O  “  Zi(i)  Zi(j)  X2(l,k)  X2(l,t)dx 

Jo 

-  -4"  WM.k.O  *  I5(i.j,k,t,4)  -  i[2I3(2,i,j  ,k,  i) 

-  I3(4,l,j,k,01  -  j  [2I3(2,  j  ,i,k,  t)  -  I3(4,j  ,i,k,  t)  ] 
+  ij [3I2(i, j ,k, t)  -  4I3(i, j ,k, t ,2)  +  I3(i, j ,k,  t ,4) ] } 


M21(i,j,k,i)  -  Zi(i)  Zi(J)  Xi(l,k)  X2(l, t)dx 

Jo 

-  _I_  {J3(k,i,j,0  -  J4(k,i,j,t,4)  -  i[2J2(2,  i,k,  j  ,  i) 

8 

-  J2(4,  i,k,  j  ,  t)  ]  -  jl2J2(2,j,k,i,0  -  J2(4,  j  ,k,  i,  i)  ] 
+  ij [ 3 J ! ( i , j ,k, O  -  4J2(i,j ,k, i ,2)  +  J2(i, j ,k, t ,4) ] } 


M22(i,j,k,i)  -  Z!(i)  Zx(j)  Xid.k)  X1(l,i)dx 

Jo 

{i2(k,i,i,j)  -  I3(k,  t  ,i. j  A)  *  i{2l1(2,i.k,  t  ,j) 

O 

-  I1(4,i,k,1>j)]  -  j[2I1(2,j,k,*,i)  *  Ii(4,j,k,t,i)] 
+  ij [3I0(i, j ,k, i)  -  4I1(i,j,k,t,2)  +  Ii(i,j,k,»,4)]J 


Ti(i,j,k) 


» 

Z!(i)  Xid.j)  X1(lfk)dx 

.  0 


-  I2(l,l. j  .k, i,l)  -  i  I0a,l,l,i,j,k) 

T2(i,j,k)  -  Z2(i)  Xjd.j)  Xi(l,k)dx 

Jo 

-  Jid.ldj.k.l)  +  i  Ji(l,l,l,j,k,i) 

T9(i, j ,k)  -  Zi(i)  X2(j,l)  X2(k,l)dx 

Jo 

-  I4(j,k,i,  1,1,1)  -  iJ2d,i,j,k,l,l) 
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-i; 


TioU.j.k)  -  Z2(i)  x2(j,l)  x2(k,l)dx 

10 

-  J3(i,j,k, 1,1,1)  +  i  J3(l.j,k,i,l,l) 
Ti4(i,j,k)  -  f  Z2(i)  X2(l,j)  X3(l,k)dx 


i 

.  0 


-  I4(l,i,j,k,l,l)  +  i  l4(l,l.i,j,k,l) 

1 


» 

.  C 


TieUJ.k)  -  I  Z!(i)  X2(1J)  X3(l,k)dx 

lo 


* 

.  I 


J5(l,i,j ,k,l,l)  -  i  J3(l,l,i,j ,k,l) 
1 


Ql(i,j,k)  -  Z3(i)  Z2(j)  Xxd.kjdx 

10 

-  Cid.i.k.j)  +  itCid.l.k.j)-  Ixd.i.j.k.l)]  -  ijli(l,l,i,k,j) 


where  Z3(i)  -  X2(l,i)  +  i[X2(i,l)  -  sin  iwx] 


,1 


Q2(i.j.k)  - 


Z3d)  Z2(j)  X2(l,k)dx 


C3(l,i,k,  j)  +  i[C3(i,l|k, j)  -  J2(l, i, j  ,k,l) ]  -  ijJ2d,l,i,  j  ,k) 


Q3d,j.k)  - 


Z3(i)  Zx(j)  Xxd.kjdx 


C2(l,i,k,j)  +  i[C2(i, l.k, j )  -  J2(l,i,k,  j  ,1)  ]  +  ijJ0dd,i,j  ,k) 
1 


Q4(i.j .k) 


Z3(i)  Zx(j)  X2(l ,k)dx 


0 


-  C4{l,k,i»  j)  +  i[C4(i,l,k,  j)  -  I3(l ,  i,  j  ,k,l)  ]  +  ijlid.l.i,  j  ,k) 
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-1 


Q5d,j,k>  -  l  z4<i>  Z2(j)  Xi(l,k)dx 

-  T2(j,i,k)  -  i[C3(k,l,i, j)  -  J2(l,j,k,i,l)]  + 

where  Z4(i)  -  Xxd.l)  -  i[Xj(l,i)  -  cos  i«x] 

.1 


Q6(iJ.k)  “ 


Z4d>  Z2(j)  X2(l,k)dx 


-  Cxd.k.iJ)  -  i[T14(j,k.i)  -  I3(l,j  ,  i .  k ,  1)  ]  +  ijl3(l,l,i,j,k) 

.1 

Qyd.j.k)  -  |  Z4(i)  Zx(j)  Xxd.kjdx 


» 

.  0 


TxCj.i.k)  -  i[C4(k,l.i.j)  -  I3d.k,i.j,l)]  -  ijlxd.l.j  ,k, i) 
_i 

QsdJ.k)  “  |  z4d>  zlU)  X2(l,k)dx 


» 

.  0 


-  C2d,k,i,j)  -  i [T16 ( j  ,k, i)  *  J4d,i.j,k,l)]  -  ijJ2d,l,j,i,k) 
1 


* 

.  C 


Ul(i,j,k)  -  I  cos  iwx  Zx(j)  Zx(k)dx 

10 

-  I5.(i,j,k,l,l)  -  jl3(l,j,i,k,l)  -  kl3(l,k,i,  j  ,1)  +  jklxdd.j  ,k,i) 

.1 

U2d»j*k)  -  |  sin  tnx  Zx(j)  Zx(k)dx 


» 

.  0 


-  J4d,j,k,l,l)  -  jJ2d,i,j.k,l)  -  kJ2(l,i,k,j,i)  ■»  jkJod.l.i.j.k) 

.1 

U3(i, j ,k)  -  j  cos  iwx  Zx(j)  Z2(k)dx 


.  0 


J4(k,i> j ,1,1)  -  jJ2(l,j,k,i,l)  +  kJ4(l,i,j,k,l)  -  jkJ2(l,l,j d,k) 
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U4(i,j.k)  -  sin  iwx  Zi(j)  Z2(k)dx 

Jo 

-  I3(i,k,j,l,l)  -  jli(l,i,j,k,l)  +  kl3(l,i,j ,k,l)  -  jkli(l,l,i,j ,k) 

U5(i,j,k)  -  I  cos  iffx  Z2O)  Z2(k)dx 

Jo 

-  I3(j ,k,i,l,l)  +  jl3(l,k,i,j,l)  +  kl3(l,j,i,k,l)  +  jkl3(l,l,i, j ,k) 

U6(i,j,k)  -  1  sin  i*rx  Z2<j)  Z2<k)dx 

Jo 

-  J2(i,j ,k,l,l)  +  jJ2(l,i.k,j ,D  +  kJ2(l,i.j ,k,l)  +  jkJ2(l.l.i*j .k) 


t 


V4(i,j,k)  -  I  X2(l,i)  Zi(j)  Zx(k)dx 

10 

J5(l,i, j ,k,l,l)  -  jJ3(l,l,j,i.k,l)  -  kJ3(l,l,k,i,j ,1)  +  jkJi(l,l,l, j ,k,i) 

.1 


V5(i,j,k)  - 


X2(l,i)  Z!(j)  Z2(k)dx 


Jo 


■I4(l.k,i,j  ,1,1)  -  Jlad.l.Jtk*1*1)  +  kl4(l,l,i,j,k,l)  -  jkl2(l,l,l.j.i.k) 
A 


V6(i,j,k) 


X2(l,i)  Z2(j)  Z2(k)dx 


J3(l. j ,k,i,l,l)  +  jJ3(l.l.k.i,j,l)  +  kJ3(l,l,j.i,k,l)  +  jkJ3(l,l,l,i, j ,k) 
1 


V7(i,j,k)  - 


X2(i,l)  Z2(j)  Z2(k)dx 


J3(i,j,k, 1,1,1)  +  jJ3(l,itk,j,l,l)  +  kJ3(l,i,j,k,l,l)  +  jkJ3(l,l,i, j ,k,l) 
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V8(i,  j  ,k)  -  X2(i,l)  ZiQ)  Z2(k)dx 

JO 

-  I4<i,k,j, 1,1.1)  *  jl2(1.1.3.k,l.l)  +  kl4(l,i,j,k,l,l)  -  jkl2(l,l. i, j ,k,l) 

V9(i,j,k)  -  X2(i,l)  ZiCj)  Zi(k)dx 

Jo 

-  J5(i,j,k,  1,1,1)  -  jj3(l.i.j,k,l,l)  *  kJ3(l,i,k, j ,1,1)  +  jkJi(l,l.i,j.k,l) 

Wjd.j.k)  -  Zi(i)  Zx(j)  Zi(k)dx 

Jo 

-  I6(i,j,k, 1,1,1)  -  il4(l,i,j,k,l,l)  -  jl4(l, j ,i,k,l,l) 

-  kl4(l,k,i, j ,1,1)  +  ijl2(l,l,i, j ,k,l)  +  ikl2(l,l,i,k,j ,1) 

+  jkl2(l,l,j,k,i,l)  -  ijkl0(l,l,l,i,j,k) 

W2(i, j ,k)  -  Z!<i)  Z2(j)  Z2(k)dx 

Jo 

-  I4(j,k,i, 1,1,1)  -  il2(l,i,j ,k,l,l)  +  jl4(l,k,i,j ,1,1) 

+  kl4(l,j,i,k,l,l)  -  ijl2(l,l,i,k, j ,1)  -  ikl2(l,l,i, j ,k,l) 

+  jkl4(l,l,i,j,k,l)  -  ijkl2(l,l,l,i, j »k) 


W3(i,j,k)  -  Zi(i)  Zx<j)  Z2(k)dx 

Jo 

-  J5(k,i,j, 1,1,1)  -  iJ3(l.i»k, j ,1,1)  -  jJ3(l,j ,k,i,l,l) 

+  kJ5(l,i,j,k,l,l)  +  ijJi(l.l,i,j,k,l)  -  ikJ3(l,l.i,j ,k,l) 
-  jkJ3(l,l,j,i.k,l)  +  ijkJi(l,l,l,i.j.k) 
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A 


W4(i,j,k) 


Pld.j.k) 


P2(iJ.k) 


P3(i.j.k) 


P4(i,j,k) 


P5(i,j ,k) 


P6<i.j .k) 


Z2(i)  Z2(j)  Z2(k)dx 

JO 

-  J3(i,j,k, 1,1,1)  +  iJ3(l,j,kti,l,l)  +  jJ3(l»i»k, j ,1,1) 

+  kJ3(l,i, j ,k,l,l)  +  ijJ3(l,l,k,i, j , 1)  +  ikJ3(l,l,j ,i,k,l) 
+  jkJ3(l,l,i, j ,k,l)  +  ijkJ3(l,l,l,i, j ,k) 

[sin  in*  -  X2(i,l)}  Z2(j)  Z2(k)dx 

Jo 

-  U6(i,j,k)  -  V7(i, j ,k) 

[sin  iirx  -  X2(i,l) ]  Z!(j)  Z2(k)dx 

Jo 

-  U4(i,j,k)  -  V8(i,j,k) 

-  [sin  iirx  -  X2(i,l)l  Zi(j)  Zx(k)dx 

Jo 

-  U2(i,j,k)  -  Vgd.j.k) 

[cos  iwx  -  X3(i,l)]  Z2(j )  Z2(k)dx 

Jo 

-  U5(i,j,k)  -  Bi(l,i,j,k) 

fl 

-  [cos  ixx  -  X3(i,l) 3  Zx(j)  Z2(k)dx 
.  0 

-  U3(i, j ,k)  -  B2(l, i, j ,k) 

-  [cos  iwx  -  X3(i,l)]  Zx(j)  Zi(k)dx 

Jo 

-  U!(i,j,k)  -  B3(l, i, j ,k) 
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The  energy  integrals  are  evaluated  in  closed  form  and  they  are  expressed 
in  terms  of  the  integrals  defined  above.  These  integrals  are  given  below. 


*  -  f  f 


f5(x#y)dxdy  -  J(l)  d(l)  - 


f5,l(x.y)f5,l(x.y)^dy  Jr2l2<l.l)l0<1^)  **  -f- 


f1  f1 

G55  -  I  fs.iO^yJfs^frty)*^  “  *2Ji(J-.i)Ji<i.i)  -  0 
Jo  Jo 

g55  —  I  j  f5,2(x>y)f5,_2<x»y)dxdy  “  *2Xo(l,l)l2(3-.»l)  “  —L 

Jo  Jo 


\ 

Jlnm  -  fi(l,y;n,m)dy  -  sin  it  [cos  n»r  I(m)  +  sin  nw  d(m)]  -  0 

Jo 


J2nm  “  f2(l,y;n,m)dy  -  cos  n*  (1-cos  7r)I(m) 


dlnm  ”  I  fl,l<X,0;n, 

Jo 


+  sin  wr  (1-sin  w)J(m)  =0  jn  **  0 

m)dx  -  jt{  [  l2(l»n)  -  nlQ(l,n)]  cos  (0) 

+  [J!(n,l)  +  nJi(l,n)]  sin  (0)}  -  0 


Kllnmpq  -  fltl(X,0;n,m)fiti(X,0;p,q)dx 
°  -  jr2[l4(n,p,l,l) 

+  npIo(l,lin,p) ] 


-  nl2(l,n,p,l)  -  pl2(l,P.n,D 


K22nmpq  -  f2,2<1»y;n’m>f2,2<1.y;P><l>dy  “  ^q*2(-l)n+Pl0(m,q) 

Jo 
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K22runpq  -  f2>3(X,0;n,m)f2(3(X,0;p,q)dx 

JO 


44 

Kllnmpq 


fl 

fl, 

Jo 


-  ?r4[  (l+n2)(l+p2)  14(11, p, 1,1)  -  (n2+p2+2n2p2)l3(n,p,l) 
+n2p2I2(n,p)  -  2n(l+p2)I2(l,n,p,l)  +  2np2Ii(l,n,p) 

-  2p(l+n2)I2(l,p,n, 1)  +  2n2pl^(l,p,n)  +  4pplQ(l,l,n,p) ] 

4(l,y;n,m)fi)4(l,y;p,q)dy  -  0 


r1  r1 

0 

F3nm  "  f3(x,y;n,m)dxdy  -  J(n)  J(m) 

•  0  «  0 


1  1 


»  r 

.  0  Jo 


-1  -1 


f4(x,y;n,m)dxdy  -  Io(l,n)Ji(l,m)  -  Ji(l,n)Io(l,m) 


*  <* 

1 

Flnm  “  f  1  f  1  (x , y ; n ,m) dxdy  -  jr{[I2(n,l)  -  nlo(l.n)]  I(m) 

Jo  Jo 


+  [JxCn.l)  +  nJi(l.n)]  J(m)> 


-  jr[J]_(n,  1)  +  nJ]_(l,n)]  J(ra)  for  m  *  0 


r1  r1 

i 

*  0  *0 


fl(2(x,y;n,m)dxdy  -  nwr[Ji(l,n)I(m)  -  nJi(l,n)J(m)] 


.1  .1 


f2ji(x,y;n, 

.  0 


0  for  n,m  *  0 


-  nwrJi(l,n)J(m)  for  m  *  0 


m)dxdy  -  x{[nJi(n,l)  +  Ji(l,n)-nJ(n) ]  I(m) 


+  [nl(n)  -  nl2(n,l)  +  Io(l,n)]  J(®)) 


r1  r1 

2 

F2nm  “  f2i2(x,y;n,m)dxdy  -  m7r{[J(n)  -  Ji(n,l)]  I(ra) 

°  °  -  [I (n)  -  I2(l,n) ]  J(m)) 


nwrl2(l.n)  J(m) 


for  n,m  t*  0 


121 


53  [  [ 

G45nm  “  f4I5<x*y‘»n*m>f5,3<x.y>dxdy 

Jo  Jo 

-  -*4{[I(l,n,l)  +  nI1(l,l,n)][J2(l,m,l)  *  m  J0(l,l,a)] 

-  [J2(l,n,l)  -  n  Jo(l*l»n)] [Ii(l,m,l)  +  mli(l,l,ni)] } 


54  f  f  53 

G45nm  “  f4t5<x,y;n,m)  f5t4(x,y)dxdy  -  G45nm 

Jo  Jo 

r1  f1 

G45nm  -  f4>5(x,y;n,m)  f5>5(x,y)dxdy 

Jo  Jo 

-  -jr4UJ2(n,l,l)  +  nJ2(l,n,l)][l3(m,l,l)  -  mlld.m.l)  ] 

-  [l3(n,l,D  -  nIi(l)n>l)][J2(m,l,l)  +  mJ2d,mtl)]> 


11 

G35nm 


■J.  1 

-f 


.2/* 

0 

1  -1 


•  I 

12 
G35nm  " 

,  0  Ji 


21 

G35nm 


22 

G35nm 


0 

a 

tt 


f3>l(x,y;n,jn)  f5(i(x,y)dxdy  -  ror2  I2(n,l)Iod»m) 

n  -  m  -  1 
all  other  n,m 

{3>l(x,r,n,m)  f5>2(x,y)dxdy  -  iwr2  Ji(l,n)Ji(m,l) 


12 

f3,2(x.y;n.“)  f5,l(x*y)dxdy  “  G35mn 


11 

f3)2(x,y;n,ni)  f5)2(x,y)dxdy  -  G35mn 
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f1  f1 

G35nm  -  f3 f 3(x,y;n,m)  f5 ( 3(x,y)dxdy  -  n2*4  I0(l,n)I0(l,m) 

Jo  Jo 


n  -  m  -  1 


all  other  n,m 


m)  f5,4(x.y)dxdy  "  G35nm 


_1  ,1 

34  f  f  33 

G35nm  “  f3i3(x,y;n,m)  f5 (4(x,y)dxdy  -  G35ran 

Jo  Jo 

G35nm  “  f3>3(x,y;n,m)  f5i5(x,y)dxdy  -  -  I2ff4Ji(n,l)Jl(m,l) 

Jo  Jo 

43  f1  f1  33 

G35nm  “  f3 ,4<x.y In»m)  f5 , 3<x.y)dxdy  “  G35mn 

Jo  Jo 

44  f1  f1  33 

G35nm  “  *3 ,4(x.y;n»m>  f5 ,4<x*y>dxdy  “  G35mn 

Jo  Jo 

45  f1  f1  35 

G35nm  “  f3 ,4<x>y ;n>m>  “  G35mn 

Jo  Jo 


44 

G35nm 


45 

G35nm 


m)  f5,4<x,y)dxdy  -  G35mn 


m)  fs.sCx-y)1131^  “  G35mn 


G35nm  -  f3 , 5(x,y;n,m)  f5i3(x,y)dxdy  -  -  nnwr4Ji(l,n)Ji(l,m) 

Jo  Jo 


54 

G35nm 


r1  f1 

f3,5(x>yin» 

Jo  Jo 


m)  f5i4(x,y)dxdy  -  G35nm 


f1  -1 

G35ran  “  f3 , 5(x,y;n,m)  f5 , 5(x,y)dxdy  -  nm7r4I2(n, l)l2(n>. D 

Jo  Jo 
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11 

G35nm  " 


12 

G45nm  ” 


21 

G45nm  “ 


22 

G45nm  “ 


33 

G45nm  ” 


34 

G45nm  “ 


35 

G45nm  “ 


f4,l<x.y;n.m)  f5,i(x,y)dxdy 


jr2([J2(n,l,l)  +  nJ2(l,n,l)]  IlU.l.m) 

-  [13(11, 1,1)  -  nlid.nil)]  Jod.l.m)} 


•  A 

f4,i<x.y:n.m)  fsdCx.y)*131^ 

0  Jo 

^{[Iid.n.l)  +  nl1(l,l,n) ]  J2(l,m,l) 

-  [J2(l,n,  1)  -  nJ0(l,l,n)  ]  lid.m,!)} 


12 

f4,2(X.y;n.“)  -5,l(x»y)dxdy  -  -  G45rtin 

0  Jo 


■  |» 

11 

f4,2(x.y;n,m)  f5>2(x,y)dxdy  -  -  G45mn 

0  .0 


* 

f4,3<x>y;n.®)  f5,3(x,y)dxdy 

0  Jo 

-  jt^{  [2nJ2(l,n,l)  -  (l+n2)J0d.l,n)l  Il(l.l.m) 
+  (2nJ1(l,n,l)  +  (l+n2)Jl(l,l,n)]  J0(l,l,m)} 


.1 

0 

.1 

0 


f4,3(x.y;n>®) 


f4,3(x.y;n,m) 


33 

f5,4(x.y)dxdy  -  G45nm 


fs.sd.y^dy 


ff4{[2nl3(n,l,l)  -  (l+n2)Ii<l,n,l) ]  J2a.ni,l) 

+  [2nJ2(n,l,l)  +  (l+n2)J2(l,n,l)  ]  Ii(l,m,l)} 
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43 

G45nm  “ 


44 

G45m 


45 

G45nm 
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